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ABSTRACT 

We present absorption line-strength maps for a sample of 18 Sb-Sd galaxies observed using 
^ \ the integral-field spectrograph SAURON operating at the William Herschel Telescope on La 

Palma, as part of a project devoted to the investigation of the kinematics and stellar popula- 
tions of late-type spirals, a relatively unexplored field. The SAURON spectral range allows the 
measurement of the Lick/IDS indices H/3, Fe5015 and Mgb, which can be used to estimate 
the stellar population parameters. We present here the two-dimensional line-strength maps for 
each galaxy. From the maps, we learn that late-type spiral galaxies tend to have high H/3 and 
low Fe5015 and Mgb values, and that the Hp index has often a positive gradient over the field, 
while the metal indices peak in the central region. 

We investigate the relations between the central line-strength indices and their correlations 
with morphological type and central velocity dispersion, and compare the observed behaviour 

■ with that for ellipticals, lenticulars and early-type spirals from the SAURON survey. We find 
r\ , that our galaxies lie below the Mg - a relation determined for elliptical galaxies and that the 

■ indices show a clear trend with morphological type. 

From the line-strength maps we calculate age, metallicity and abundance ratio maps via a 
comparison with model predictions; we discuss the results from a one-SSP (Single Stellar 
Population) approach and from a two-SSP approach, considering the galaxy as a superposi- 
tion of an old (« 13 Gyr) and a younger (age ^ 5 Gyr) population. We confirm that late-type 
galaxies are generally younger and more metal poor than ellipticals and have abundance ra- 
tios closer to solar values. We also explore a continuous star formation scenario, and try to 
recover the star formation history using the evolutionary models of Bruzual & Chariot (2003), 
assuming constant or exponentially declining star formation rate (SFR). In this last case, fix- 
ing the galaxy age to 10 Gyr, we find a correlation between the e-folding time-scale r of the 
starburst and the central velocity dispersion, in the sense that more massive galaxies tend to 
have shorter r, suggesting that the star formation happened long ago and has now basically 
ended, while for smaller objects with larger values of r it is still active now. 

Key words: galaxies: bulges - galaxies: evolution - galaxies: formation - galaxies: kinematics 
and dynamics - galaxies: spiral -galaxy: structure 
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1 INTRODUCTION 

Over the last few decades, stellar population synthesis has be- 
come one of the most popular and powerful techniques to study 
the formation and evolution histories of galaxies too distant to be 
resolved into individual stars. In particular, the measurement of 
absorption line-strengths combined with stellar population mod- 
els has been used to investigate the luminosity-weighted age, 
metallicity and abundance ratios in integrated stellar populations. 
Extensive work has been carried out on stellar populations in 
early-type galaxies, mainly using long-slit spectra (see for exam- 
ple [Carollo, Danziger, & Buson 1993, Davies, Sadler & Peletier 
1993, Fisher, Franx & Illineworth 1995, Fisher, Franx & lUing- 



worth 1996, Longhetti et 



in eworth 1993, Fisher, Franx & lltmg- 
^ 1 19981 |j0rgensenlll999l. iKuntschneil 



bOOd. ISanchez- Blazquez et al.H2007l) and recently also integral 
field spectroscopy (Kuntschner et al. 2006, hereafter Paper VI, Mc- 
Dermid et al. 2006, hereafter Paper VIII). The main advantage 
brought by the integral-field observations is that one can obtain 
two-dimensional maps and easily identify extended structures. In 
particular, the two-dimensional coverage gives the possibility to 
straightforwardly connect the stellar populations with the kinemat- 
ical structures and therefore represents a powerful tool towards a 
comprehensive understanding of the structure of the galaxy. 

Some literature is also available on absorption line- 
strengths in early-type spiral galaxies, mainly SO and Sa galax- 
ies; (see for example Jablonka, Martin & Arimoto 1996, 
Idiart. de Freitas Pacheco & Costj[l993. IProctor & Sansomll2002l 



Afanasiev & Sil'chenkdl2005l. and Peletier et al. 2007, hereafter 
Paper XI). In particular. Paper Xj presents two-dimensional ab- 
sorption line-strength maps of 24 spiral galaxies, mostly of type 
Sa, observed with the integral-field spectrograph SAURON (Bacon 
et al. 2001, hereafter Paper I); these reveal young populations in 
Sa galaxies, possibly formed in mini-starbursts in structures such 
as nuclear discs, circular star forming rings and bars. These mini- 
starbursts generate a large scatter in index - index diagrams, larger 
than observed for elliptical galaxies. As a consequence, there is 
a large range in luminosity-weighted ages. Different star forma- 
tion modes -starbursts or quiescent- are also reflected in the wide 
range observed in the abundance ratios. All the galaxies of Paper 
XI lie on or below the Mgfo - a relation for ellipticals determined 
by l0rgensen, Franx & Kjaergaard (1996). The authors argue that, if 
one considers that relation as valid for old galaxies, then the early- 
type spirals present a scatter in age, with the biggest scatter for the 
smallest velocity dispersions. 

Spiral galaxies towards the end of the Hubble sequence are, 
on the other hand, still poorly-studied objects, due to the difficulty 
of obtaining reliable measurements in these low-surface bright- 
ness objects, full of dust and star forming regions, whose spec- 
tra can be dramatically contaminated by emission-lines. In par- 
ticular, the study of their stellar populations is extremely compli- 
cated. An approach based on broad-band colours is basically in- 
validated by the huge amount of dust, which is very difficult to 
take into account, especially when working with data from ground- 
based telescopes. A spectroscopic approach, relying for example on 
stellar absorption-strength indices, has instead to deal with the al- 
most ubiquitous presence of gas with consequent emission, which 
makes it necessary to accurately remove the emission-lines from 
the spectra before performing any population analysis. Only in 
the last decade several HST-based imaging surveys have targeted 
these late-type gala xies. Particularly relevant are the papers by 
ICarollo et alj \\99%. Carollo & Stiavelli 1998, Carollo, StiaveUi 



& Mack 1998. ICarollol ( Il999l) , ICarollo et al.l ( l2002h . iBoker et al.l 



| |2002|) .I Laine et al. I( l2002l) . where the authors reveal the presence of 
a variety of structures in the inner regions, such as bulges, nuclear 
star clusters, stellar discs, small bars, double bars and star forming 
rings whose fo rmation and evolutio nary pattern are not properly 
understood yet. ICaroUo et al.l ( |2006|) also performed a population 
analysis on nine late-type spiral galaxies (types between Sa and Sc) 
on the basis of HST ACS and NICMOS optical and near-infrared 
colours. The high spatial resolution of their space-based data al- 
lowed them to mask dust features and measure the properties of the 
stellar populations from the colours; this method, together with the 
high spatial resolution provided by HST, produced results that do 
not suffer much from the effects of dust; only a smooth dust dis- 
tribution would not be detected. From their analysis, they found a 
large range in the colour properties of bulges in late-type spirals, 
reflecting a large range in their stellar population properties. From 
a comparison with population models, they concluded that in about 
half of their bulges the bulk of the stellar mass formed at early 
epochs (more than 50% of their mass was formed more than 9 Gyr 
ago), and that in the other half a non-negligible fraction of the stel- 
lar mass (up to 25%) was formed recently, in the last « 3 Gyr. 

Until now, only little spectroscopic counterpart for the quoted 
imaging was available, due to the men tioned observational difficul- 
ties. The published work ( iMatthews & Gallagher 20 ()3, Zaritsky , 
Kennicutt & Huchr a 1994, lB6ker'eraL,200L .Rossaet alj|2003l 
IWalcher et a"i]|2005h mainly refers to the emission-line properties, 
to the characteristics of the HII regions or to the nature of the in- 
nermost component. In particular, no study has addressed yet in a 
systematic way the stellar populations of late-type spiral galaxies 
with a two-dimensional coverage. Existing studies on abundances 
concern topics such as elemental abundan ces in HII region s and 
the distribution of oxygen within discs (Zarits ky at alj[l994h . not 
properties inferred from the integrated stellar absorption spectrum. 
As for the absorption-lines properties, there is ongoing work by 
MacArthur, Gonzalez & Courteau (2007), who analyse long-slit 
spectra (Gemini - CMOS) of eight spiral galaxies of type between 
Sa and Scd on a wide spectral range, measurin g many absorption 
indices out to 1-2 disc scale-lengths. Also, Moorthv & HoltzmanI 
( l2006h presented line-strengths in the bulges and inner discs of 
38 galaxies of morphological type between SO and Sc for which 
they acquired long-slit spectra on a broad spectral range including 
Balmer lines, Mg and Fe features. They concluded that the cen- 
tral regions of bulges span a wid e range both i n SSP metallicity 
and age, confirming the results of ICarollo et al.l ( |2006|) . They also 
found that luminosity-weighted metallicities and abundance ratios 
are sensitive to the value of the central velocity dispersion and of 
the maximum disc rotational velocity, and that red bulges (defined 
as those with B~K > 4) of all types are similar to luminous el- 
lipticals and obey the same scaling relations; for the blue bulges 
(B— K < 4) they observed instead some differences. They also ad- 
dressed the radial variations of the population parameters; in most 
cases, they measured negative metallicity gradients with increas- 
ing radius, in the bulge-dominated region; positive gradients in age 
were found only in barred galaxies. 

As mentioned, there is still a lack of work on populations 
in late-type spirals on a full two-dimensional field: we therefore 
started a project on a sample of 1 8 late-type spiral galaxies, with the 
purpose of investigating the nature of their inner regions, address- 
ing the bulge and disc formation and evolution, the interconnection 
between stellar and gaseous components and the star formation his- 
tory. Given the high complexity in the inner regions of these ob- 
jects, integral-field spectroscopy has to be preferred with respect to 
long-slit, providing a full picture of the kinematics and populations 
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on a two-dimensional field of view. We observed our 1 8 late-type 
spiral galaxies using SAURON, an integral-field spectrograph that 
was built for a representative census of ellipticals, lenticulars and 
early-type spiral galaxies, to which in the rest of this Paper we will 
refer as 'the SAURON survey'. Our project can be regarded as an 
extension of the SAURON survey towards later-type objects. First 
results based on our data are presented bv lOanda et alj ( l2006f) and 
focus on the stellar and gas kinematical maps. The main findings 
of that Paper are that in many cases the stellar kinematics suggests 
the presence of a cold inner region, as visible from a central drop in 
the velocity dispersion, detected in about one third of the sample; 
that the ionized gas is almost ubiquitous and often presents more 
irregularities in its kinematics than the stellar component; and that 
the line ratio [OIIIJ/H/3 assumes often low values over most of the 
field, possible indication for wide-spread star formation, in contrast 
to early-type galaxies. 

As a following step in our project, we subtracted the 
e mission-lines f rom o ur SAURON spectra, measured as explained 
m iGanda et alj ( l2006h , and on the emission-cleaned spectra de- 
rived absorption line- strength indices in the Lick/IDS system 



orpt] 

jpaber et alJll985LlGorgas et alJI 19931 . IWorthev et al.lll994l ). which 
allows an easier comparison with existing data and with other 
classes of objects. Here we present and analyse the absorp- 
tion line-strength maps. The data and maps presented in this 
Paper will be made available via the SAURON WEB page 



[http : / /www .strw.leidenuniv.nl/sauron/ 

The Paper is structured as follows. Section 2 briefly describes 
the sample selection, observations and data reduction. Section 3 
reviews the analysis methods. In Section 4 we present the line- 
strength maps for all the galaxies in our sample. In Section 5 we 
investigate the correlations between different indices and between 
indices and other galactic parameters, such as the morphological 
type and the central velocity dispersion, mainly focussing on the 
central values. In Section 6 we present our estimates for the ages 
and metallicities, obtained by comparison of our observed line- 
strength indices with model predictions; we also explore the star 
formation history and give estimates for the time-scale r, in an ex- 
ponentially declining star formation rate scenario. Section 7 sum- 
marises the main results. In Appendix A we deal with the radial 
variations of the indices and the stellar population parameters, pre- 
senting azimuthally averaged radial profiles. Appendix B shortly 
describes the individual galaxies and their peculiarities, giving a 
brief summary of relevant literature information. For more litera- 
ture reviews on these objects, we suggest the reader to check Sec- 
tion 6.4 in .Ganda et al., C2006) . 



2 SAMPLE SELECTION, OBSERVATIONS AND DATA 
REDUCTION 

The galaxies wer e optically selected (Bt < 12.5, according to 
the values given in lde Vaucouleurs et al.l 19911 hereafter RC3) with 
HST imaging available from WFPC2 and/or NICMOS. Their mor- 
phological type ranges between Sb and Sd, following the classifi- 
cation given in NEE0 (from the RC3). Galaxies in close interac- 
tion and Seyferts were discarded. The resulting sample contains 18 
nearby galaxies, whose main properties are listed and illustrated re- 
spectively in Table 1 and Fig. 1 in Ganda et al. (2006). The NGC 
numbers and some characteristics are reported also in the tables in 
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the present Paper. Observations of the 18 sample galaxies were car- 
ried out during 6 nights in January 2004, using the integral-field 
spectrograph SAURON attach ed to the 4.2 -m WHT. The exposure 
details are listed in Table 2 in lGandaetalj f2006). 

We used the low spatial resolution mode of SAURON, giv- 
ing a field-of view (FoV) of 33" x 41". The spatial sampling of 
individual exposures is determined by an array of 0'.'94 x 0'.'94 
square lenses. This produces 1431 spectra per pointing over the 
SAURON FoV; another 146 lenses sample a region 1^9 away from 
the main field in order to measure simultaneously the sky back- 
ground. SAURON delivers a spectral resolution of 4.2 A FWHM, 
and its maximum spectral coverage is 4800-5380 A; the wavelength 
interval common to all lenses is the narrow spectral range 4825- 
5275 A (1.1 A per pixel). This wavelength range includes a num- 
ber of important stellar absorption lines (e.g. H/3, Fe5015, Mgb), 
and potential emission-lines as well (H/3, [Q UI], [NI ]). For a more 
exhaustive description of the instrument, see lPaper"^ 

Data reduction was carried out using the dedicated software 
XSAURON developed at Centre de Recherche Astronomique de 
L yon (CR AL) . For details on the d ata reduction we refer the reader 
to Paper 1 and lGanda etakl ( |2006|) . For an explanation of the spec- 
trophotometric calibration we recommend the reader to study Paper 
VI, in particular its Section 3.1. 



3 ANALYSIS AND METHODS 

The fully reduced and calibr ated spectra were first spat ially binned 
using a Voronoi tessellation dCappellari & CoDinll2003t) in order to 
reach a minimum signal-to-noise ratio of « 60 per spectral reso- 
lution element per bin. The kinematics of stars and gas were then 
extracted from the bi nned spectra via the pixe l -fittin g methods ex- 
tensively described in lCappellari & EmsellemI ( l2004h . Emsellem et 
al. (2004, hereafter Paper III), and Sarzi et al. (2006, hereafter Pa- 
per V), t hat we already app lied to calculate the kinematical maps 
shown in lGanda et alj(l200q). The spectra were fitted using the SSP 
models from VazdekisI ( Il999l) . We selected a library of 48 models 
evenly sampling a wide range in both age and metallicity (0.50 ^ 
Age s; 17.78 Gyr, -1.68 ^ [Fe/H] < +0.20). We used this SSP h- 
brary and additional Gaussian templates reproducing the emission- 
lines to fit the observed spectra, separating the emission -lines from 
the underlying absorption line spectrum (see IPaper Vl for details 
on the procedure). We notice that for our objects this separation is 
crucial if we want to measure absorption strengths, since in spiral 
galaxies the gas is almost ubiquitous and emission in several cases 
dominates the spectrum, filling in completely the absorption fea- 
tures, as shown for example in the cases of some central aperture 
spectra in the top row of Fig. [T] (black lines); in the same Figure we 
overplot the best-fit to the spectra obtained with our method (red 
lines) and the emission-removed spectra (blue lines), illustrating 
the reliability of the applied gas-cleaning procedure. The bottom 
row of Fig.[T]presents similar plots for bins located « 10" from the 
galaxy's centre, showing that our fits and emission correction are 
reliable also outside the central region, given the fact that a mini- 
mum signal-to-noise of 60 is garanteed by the spat ial binnin g. The 
re ader can find more examples of this in Fig. 1 in IPaper Vl Fig. 3 
in iGanda et alj ( l2006h . and Fig. 3 in Falcon-Barroso et al. (2006), 
hereafter Paper VII. 

We also notice that in two cases (namel y, NGC 2964 and 
NGC 4102) we detected regions of activity (see lcanda et al.ll2006l 
for a description), where the spectra are characterised by double- 
peaked emission-lines due to the presence of multiple velocity 
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Figure 1. Top row: central aperture spectra of some representative galaxies in our sample: from left to right, we plot NGC 1042, 2964 and 4487. Along the 
vertical axis, the flux is in units of 10"^'' erg cm~^ s~^. The black lines -hardly distinguishable from the other ones- represent the observed spectra; the red 
ones the best-fit and the blue ones the spectra after emission removal. Bottom row: as in the top row, but for bins located ~ 10" away from the galaxy's centre. 



components in the gas, because of the presence of a jet. In those 
regions our method of fitting the fines with single Gaussians fails 
to reproduce accurately the line profile, but since this is limited to 
only a few spectra, we did not modify the method. In any case, 
this does not have any impact on our following analysis, since the 
double-peaked emission-lines do not contaminate the spectra defin- 
ing the central aperture (see Section|5]l on which most of this study 
is focused on. 



3.1 Line-strengths measurement and FeSOlS correction 

On the emission-removed spectra we calculated line-strength in- 
dices on the Lick/IDS system. For det ails on th e actual compu- 
tation of the line-strengths, we refer to IPaper v3 . The calibration 
to the Lick/IDS system takes into account several effects: the dif- 
ference in spectral resolution between the Lick/IDS system and 
the SAURON instrumental setup, the internal velocity broadening 
of the observed galaxies and small systematic offsets rising from 
differences in the shape of the continuum, due to the fact that the 
original Lick/IDS spectra have not been flux calibrated. In our ob- 
serving run, we obtained spectr a for 21 di fferent stars in common 
with the Lick/IDS stellar library jWorthev et al. 1994), for a total of 
39 repeated observations. Fig. |2] shows the difference between the 
Lick/IDS measurements and ours (for stars with repeated observa- 
tions, we considered the average of the different values available). 
We calculated the mean offsets to the Lick/IDS system and their 



dispersion using bi-square weighting in order to minimize the in- 
fluence of the outliers, after removing from the sample two stars of 
spectral type M6III and M7III, whose measured metal indices dif- 
fer more than 2 A from the Lick measurements. The adopted offsets 
that we applied to all our measurements are: —0.06 A (± 0.06 A) 
for H/3, -1-0.32 A (± 0.15 A) for Fe5015 and -1-0.15 A (± 0.05 A) 
for Mgb. 

From a visual inspection of our maps of the Fe5015 index, we 
saw that our measurements tend to assume high values in the outer 
parts, especially in the bottom-right corner, revealing a systematic 
effect that makes the maps look asymmetric. The most striking 
manifestation of the problem is this 'high-iron comer', but there 
could be as well regions where the index is too low. A detailed in- 
spection of the spectra showed that the cause is a wrong shape of 
the continuum in the Fe5015 spectral region, as proved by the fact 
that we are unable to fit our twilight spectra using our SSP tem- 
plate library: the Fe5015 region is not matched by the templates, 
and the Fe5015 map calculated on the twilight spectra is not flat 
over the field, showing a range of about 1.6 A from minimum to 
maximum value. This continuum problem could be due to instru- 
mental instabilities, to the high sensitivity of the grism to the inci- 
dence angle of the light, to the imperfect extraction of the spectra, 
to the misali gnment of the spec tra with respect to the columns of 
the CCD (see lOanda et al.l2006l) . We can exclude, instead, relations 
with problems in the gas removal, since the Fe5015 measurement is 
affected also in the twilight frame, where the spectra are gas-free. 
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Figure 2. Differences between the Lick/IDS and our measurements for the 21 stars in common, for the three indices. For two stars (HD18191 and HDl 14961, 
spectral type respectively M6III and M7III) our measurements of the metal indices differ dramatically from the Lick/IDS measurements; therefore we excluded 
them from the determination of the offsets and plotted them with a different symbol (open triangles). The overplotted dashed line shows the adopted offset, 
derived by a bi-weight estimator In the bottom right corner we place a typical errorbar. 
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Figure 3. Maps of the Fe5015 index for the twilight frame, before and after correction, in A. The two maps are plotted with the same cuts, indicated on their 
side, together with the colour bar. In the corrected frame, the spatial gradient is largely removed. 



In addition to this, continuum problems should not influence the 
gas-cleaning procedure: in the fitting, a polynomial is included in 
order to account for differences in the flux calibration between the 
observed spectra and the model ones. 

Unfortunately, we could not take this effect into account prop- 
erly at the data reduction stage, and decided to apply a correction to 
the affected spectra before measuring the Fe5015 index. A problem 
related to the shape of the continuum was experienced also for the 
early-type spirals of the SAURON survey and affected the H/3 mea- 
surements, due to the closeness of the spectral feature to the edge 
of the spectral range; despite the fact that the problem is different, 
we worked out a correction similar to the one that was there applied 



in the case of H/3 and described in lPaper v3 (Section 3.1.2). 

For each galaxy, we averaged the spectra within an aperture of 
10" radius, obtaining a 'global spectrum' for that galaxy. We then 
fitted this global spectrum using a linear combination of SSP tem- 
plates, masking the spectral regions possibly affected by emission 
lines. In this way we obtained for each galaxy a 'global optimal 
template' that approximates very well the general spectral features 
of that galaxy. Then we considered the two-dimensional galaxy's 
datacube, spatially binned and cleaned from emission as explained 
in the previous Sections; for each bin, we fitted the cleaned spec- 
trum using the 'global optimal template' together with an 1 1-th or- 
der multiplicative polynomial, over the whole spectral range. The 
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Figure 4. Maps of the Fe5015 index for NGC488, before and after correction, in A. Tlie wliite asterisks in the right-hand map mark the position of the 
centroids of the bins to which we have applied the correction. Overplotted on the maps are the isophotal contours. The two maps are plotted with the same 
cuts, indicated on their side, together with the colour bar; they are oriented with the horizontal and vertical axis aligned, respectively, to the long- and to the 
short- axis of the SAURON field. 



fit determined the best-fitting polynomial continuum, from which 
we then removed the hnear component (which does not change the 
line indices); we considered the residuals and calculated their RMS 
variation over the wavelength range of Fe5015. For bins covering 
four or more single-lens spectra and with an RMS greater than 0.01, 
we used the fitted polynomial as a correction to the spectral shape 
by dividing it through the emission-cleaned spectrum. For these 
bins, we finally measured the Fe5015 index on the corrected spec- 
trum. 

This correction depends critically on the assumption that the 
optimal template is an accurate representation of the true galaxy 
spectrum. Whilst this is generally true, in the central, sometimes 
bulge-dominated, regions of our galaxies, non-solar abundance ra- 
tios can create a fit residual similar in size to the systematic ef- 
fect we wish to remove. By correcting only bins containing four or 
more lenses, we avoid these central parts, since the data there are 
unbinned. Correcting large bins is also important because of their 
larger contribution to the spatial appearance of the maps. 

We verified the Fe5015 correction using exposures of the twi- 
light sky. As we mentioned above, in the twilight maps of Fe5015, 
a spatial gradient was also visible, as is apparent in the galaxy data. 
In the uncorrected twilight, the field centre gave values of Fe5015 
close ~ 4 A, consistent with a typical G2V - G4V sun-like spec- 
trum, showing that the field centre is not significantly affected by 
the problem. After applying our correction techniquqj, the gradient 
was reduced, and the global level remained around 4 A. We show 
this in Fig. (3) where we display the Fe5015 map calculated on the 
twilight frame, before and after applying the correction. Similarly, 
in Fig.|4]we show the Fe5015 map for NGC488 before and after 
correction. White asterisks mark the centroids of the bins that have 
been corrected. The correction clearly removes most of the asym- 
metry, proving the usefulness of the method applied. 

Figure|5]presents for each galaxy the map of the RMS devia- 
tions of the fitted continuum after removal of a linear slope. We plot 
in black the bins that do not meet the requirements for the correc- 
tion, for which we have artificially put the RMS to (for plotting 

^ applied in this case to all single-lens spectra with an RMS greater than 
0.01: given the high signal-to-noise of the twilight frame, no spatial binning 
is necessary 



purposes only); in these bins the Fe5015 index is calculated on the 
original emission-cleaned spectra. The maps are oriented with the 
horizontal and vertical axis aligned, respectively, to the long- and to 
the short- axis of the SAURON field. The extension of the corrected 
area is mainly related to the signal-to-noise of the data: galaxies 
with poor signal-to-noise will require higher binning and the big 
bins will undergo the correction. 

As for the uncertainties on the measured indices, the most 
worrying factor is the separation of emission and absorption lines. 
Simulations presented in Appendix A of Paper V show that the ac- 
curacy in recovering the emission-line fluxes does not depend on 
the A/N ratio between the line amplitude and the noise level in 
the stellar continuum, but only on the signal-to-noise level in the 
stellar continuum itself. These simulations were meant to investi- 
gate the accuracy of emission-line fluxes measurements for ellipti- 
cals and lenticulars. In the case of signal-to-noise = 60, the typical 
uncertainties in the line fluxes obtained in Appendix A of Paper 
V translates in errors in the equivalent width of the emission line 
of 0.0 8 A, and in similar errors in the absorption line indices. 
IPaper xj presents absorption line-strengths for early-type spirals 
observed with SAURON; in some of the galaxies studied there, the 
emission lines are stronger by a factor up to 100, relative to the ab- 
sorption lines, than in the ellipticals or lenticulars of Paper V with 
the strongest emission-lines. Therefore, for the objects in Paper XI 
the situation is much closer to what happens in our spirals. There, 
the authors report results of similar simulations of the reliability of 
the gas cleaning, exploring a range in A/N ranging up to 100. The 
results are such that the errors in the absorption line indices do not 
vary appreciably with increasing A/N. When coming to the data, 
one has to consider that the spectra are packed close together, so 
that the absorption lines are affected by neighbouring spectra; de- 
spite the reduction tries to minimize this contamination, there might 
be some residual effects, particularly in presence of strong emis- 
sion lines. Additional errors can derive from template mismatch in 
the fitting procedure and from the continuum correction described 
above. Therefore, throughout the rest of the Paper we adopt a con- 
servative value of 0.2 A for the uncertainties in H/3 and Mgi and 
0.3 A in Fe5015. The same figures were used in Paper XI. 
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Figure 5. Galaxy by galaxy maps of the RMS deviations of the fitted polynomial continuum after removal of a linear slope, in the correction procedure for the 
continuum problem affecting the Fe5015 measurement; in the bins in the black areas the connection will not be applied. All the maps are plotted with the same 
spatial scale and oriented with the horizontal and vertical axis aligned, respectively, to the long- and to the short- axis of the SAURON field. 
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4 LINE-STRENGTH MAPS OF 18 LATE-TYPE 
GALAXIES 

Figures I6al6il present the absorption line-strength maps of the 18 
galaxies, ordered by increasing NGC number. For each galaxy, in 
the first row we show the total intensity obtained by integrating 
the full SAURON spectra in the wavelength direction, the galaxy 
NGC number (and the UGC also, when available), and the two- 
dimensional map of the characteristic time-scale of star formation 
r, in Gyr (see Section [63). Overlapped in pink on the intensity 
map are the contours delimiting the outer limits of the so-called 
'bulge' and 'disc regions', defined in Section[5] In the second row 
we present the two-dimensional maps of H/3, Fe5015 and M-gb, 
where the measured indices have been calibrated to the Lick/IDS 
system and are expressed in equivalent widths and measured in A. 
The maps of Fe5015 have been corrected as explained in Section 
13.11 the corrected bins are marked with a white dot. In the third 
row, we present the maps of age (expressed in units of decimal 
logarithm and measured in Gyr), metallicity (in units of decimal 
logarithm, with the solar metallicity as zero point) and abundance 
ratios obtained with the one-SSP approach (see Section [6Tt . Over- 
plotted on each map are the isophotal contours. The maps are all 
plotted with the same spatial scale, and oriented with the horizon- 
tal and vertical axis aligned, respectively, to the long- and to the 
short- axis of the SAURO N field for presentati on purposes; the ori- 
entation is the same as in lGanda et al.l ( llOOd) and as in Figures |4] 
and|5]in this Paper. The relative directions of North and East are in- 
dicated by the arrow above the galaxy's name. The maximum and 
minimum of the plotting ranges and the colour table are given in the 
tab attached to each map. To make the comparison among different 
objects more immediate, we adopted the same plotting ranges for 
all of the galaxies, but in Appendix IbI we will enlighten the spatial 
structures within the single maps, when noticeable, by inspecting 
the maps plotted on a 'galaxy-by-galaxy range'. 
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Figure 6a. Maps for NGC 488 and NGC 628. For each galaxy; first row: reconstructed image; NGC and UGC numbers and time-scale r for star formation in 
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5 CORRELATIONS BETWEEN LINE-STRENGTH 
INDICES AND OTHER GALACTIC PARAMETERS 

In this Section we analyse the dependence between the line- 
strength indices and other galactic parameters; we mainly focus on 
central properties. For each galaxy we defined three regions: the 
central circular aperture of l'.'5 radius (corresponding to the inner- 
most nine single-lens spectra); a so-called 'bulge region', defined 
as the region within which the light due to the fitted exponential 
disc falls below 50% of the total light, with the exclusion of the 
central aperture; and a so-called 'disc region', a 3"-thick elliptical 
annulus external to the 'bulge region|f|. We averaged the spectra 
in these regions before spatial binning, and on the resulting spec- 
trum calculated the line-strength indices, after emission removafl 
The actual numbers for the definition of the bulge and disc regions 
(semi-major axis of the ellipses, position angles and ellipticities) 
were taken from our own photometric analysis and bulge-disc de- 
composition based on archival NIR images (Ganda et al., in prepa- 
ration). There we model the galaxies as an exponential disc and a 
Sersic (1968) bulge (which often turns out to be very tiny, with ef- 
fective radius smaller than ~ 5") and apply a two-dimensional de- 
composition from w hich we can estimate the bulge extension (cfr. 
lNoordermeerll200^ for the method and Ganda et al., in preparation, 
for the actual analysis of the profiles). For clarity, in the top left map 
in Figures [6a|6i| we overplot in pink the elliptical contours defining 
the outer limits of the 'bulge region' and of the 'disc region'. In any 
case, we want to warn the reader that the 'bulge' and 'disc' regions 
here introduced are mainly a tool for us to investigate radial vari- 
ations, and do not refer to distinct components, since the so-called 
'bulge' is affected by a huge contamination from the disc, and an 
accurate decomposition is beyond the scopes of the present Paper, 
and will be addressed in detail in the mentioned paper in prepara- 
tion. 

In Table[T]we list the morphological type (from the RC3) and 
the stellar velocity dispersion, with its associated error, measured 
on the central aperture spectra, since we will make use of these 
quantities in the following Sections and Figures; the errors on the 
velocity dispersion are based on the scatter of the velocity disper- 
sion values that we measure for the single-lens spectra within the 
central aperture, and do not take into account template mismatch. 
Table |2] lists the index values calculated, galaxy by galaxy, on the 
central aperture and on the bulge and disc regions. 



5.1 Index - index relations 

In order to derive 'zero-order' estimates on the luminosity weighted 
ages and metallicities of our galaxies, we show in Fig. |7] an 
age/metallicity diagnostic diagram. Here we use the abundance ra- 
tio insensitive index 



[MgFe50] = 



0.69 X Mgi + Fe5015 



(1) 



^ The inner semimajor axis of the 'disc region' coiTesponds to the outer 
semimajor axis of the 'bulge region' plus 2". A slightly different definition 
of the disc region has been adopted in the case of NGC 628, in order to avoid 
contamination from a foreground star situated ss 13" south of the centre. 
* For the aperture measurements, we never applied the continuum correc- 
tion to the Fe5015 index described in Section IxT] in several cases, even 
the disc region is affected only partially by the correction; in addition to 
this, the bulge and disc region measurements are obtained averaging spec- 
tra from spatially sepai'ated ai'eas, which helps in minimising the impact of 
the continuum problem in the Fe5015 spectral region. 
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Table 1. Hubble type (RC3 through NED); numerical morphological type 
(from RC3, indicated as 'T-type' in some of the following pictures); stellai' 
velocity dispersion (in km s~^) and its associated error, also in km s~^, 
measured on the central aperture spectrum, for all galaxies. 




1 2 3 4 5 6 
[MgFeSO] (A) 

Figure 7. H/3 index against [MgFe50] (in A), for three samples of galaxies 
observed with SAURON: the black symbols represent the central apertures 
of our late-type spirals, the red symbols are the E and SO (central aperture of 
1"335 radius) and the blue ones th e Sa galaxies ( central aperture). The solid 
lines represent SSP models from iThomas et al J 1(2003,), limited to models 
with solar abundance ratios; ages of the models vary from 0.6 to 15 Gyr 
(top to bottom); metallicity (in decimal logarithm) from —2.25 to +0.35 
(left to right). Representative error bars are added at the bottom right of the 
panel; the black one refers to both the early- and late-type spirals samples, 
while the red one, smaller, refers to the E/SO galaxies. 
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Table 2. Line-strength indices calculated on the central apertures and bulge and disc regions, defined as explained in Section|5]for the 18 galaxies. The indices 
are intended as equivalent widths and expressed in A. 
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Figure 8. Left panel: index - index diagram showing Fe5015 against yigb for the central aperture of our late-type spirals (black symbols), while the red 
symbols are the E and SO (within a 1" 335 central aperture) and the blue ones the Sa (central aperture), observed with SAURON. The solid lines represent the 
SSP models from lThomas et alj j2003l) . limited to models with solar abundance ratios; the dashed fines correspond instead to the models with [Mg/Fe]=0.5. 
Middle panel: same as the previous, for the 'bulge region' apertures of the 18 late-type spirals only. Right panel: same as the previous, for the 'disc region' 
apertures of the 18 late-type spirals only; in the middle and right panels we only plot the model grid for solar abundance ratios. All the indices are expressed 
as equivalent width and measured in A; the three panels share the same vertical axis. Representative error bars are added at the bottom right of each panel; the 
black one refers to both the early- and late-type spirals samples, while the red one refers to the E/SO galaxies. 



(for details on the index definition see Kuntschner et al., in prepa- 
ration) as metallicity indicator and H/3 as age indicator. In Fig. |7] 
we plot the indices calculated on the central I'.'S aperture of our 
galaxies (black symbols). We overplot with red symbols the ellip- 
tical and lenticular g alaxies fro m the SAURON survey (based on 
the data presented in IPaper v3) . with reference to measurements 
within apertures of radius 1'.'335, that the authors of that pa per ex- 
tracted in order to match the Lick aperture (see IPaper VI for de- 
tails); the blue symbols represent central ap ertures of l"2 radius for 
the early-type spirals from the same survey jPaper x3) . Overplotted 
in the Figure are also the models of Thomas, Maraston & Bender 
(2003), Umited to the models with [Mg/Fe]=0; they ar e labelled 
with the corresponding age and metallicity. As noted in lPaper 5c3. 



there is a smooth transition between the ellipticals/lenticulars and 
the early-type spirals, with most E/SO populating the region of the 
old, metal-rich models, and the Sa spanning a larger range in age 
and reaching lower metallicities. Our 18 later-type spirals occupy 
a region of the diagram that has little intersection with the E/SOs, 
having generally younger ages and being clearly more metal poor. 
We also notice that the parameter region occupied by our galaxies 
is completely included in the one spanned by the Sa galaxies. We 
warn the reader that the interpretation in terms of ages and metal- 
licities is here done on the basis of SSP models, which might not 
provide a good description of spiral galaxies, for which the star for- 
mation history could be characterised by several bursts spread over 
time or by a constant star formation rate. A more detailed analysis 
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of the age, metallicities and the star formation history is presented 
in Section|6] 

In Fig. [S] the left panel shows the line index Fe5015 against 
MgZ> for the central aperture of our 18 galaxies (black symbols); 
as in Fig. [T] we show also the points corresponding to the E/SO 
(red symbols) and Sa (blue) of the S AURON survey. The solid 
lines represent the SSP model grid from lThomas et all ^20031) with 
[Mg/Fe]=0; the dashed lines correspond i nstead to the Thomas 
models with [Mg/Fe]=0.5. As also shown in lPaper 5c3 with a simi- 
lar figure, it is clear that most of the ellipticals can be fitted better by 
models over-enhanced in a-elements, while the Sa bulges present a 
range in enhancement. Our later-type spirals instead are in all cases 
compatible with models with solar-abundance ratios. In the middle 
and left panels we show, only for our sample of late-type spirals, the 
measurements for the previously defined bulge and disc regions re- 
spectively in the (Mgii, Fe5015) plane, along with the models from 
[Thomas et al. ( 2003), drawn only in the case of solar abundance ra- 
tios. We see that our spirals are compatible with solar abundance 
ratios, within the uncertainties; we do not detect relevant global 
variations in the appearance of the diagram when considering the 
different regions (central aperture, bulge and disc region): the sin- 
gle galaxies move in the diagram, but they remain in the region 
of the models with solar abundance ratios. This differs from what 
seen in the early-type spirals of the SAURON survey: there, away 
from the centre (~ 10" from the centre) the galaxies have generally 
higher [Mg/Fe] than in the centre. Super-solar abundance ratios are 
possibly caused by the presence of populations with several ages. 
An increase of [ Mg/Fe] when moving outwards is det ected also by 
Ijablonkj ( |2007|) and by iMoorthv & HoltzmanI ( l200d ) for samples 
of galaxies of type varying between SO and Sc. This might be due 
to a change from SSP-like stellar populations in the cent re to pop- 
ulation s with a large range in age outside. In particular, Ijablonkal 
( l2007h measures small increments in abundance ratio from the cen- 
tral region to the bulge effective radius; th ese gradients stay rather 
constant among bulges; lMoorthv & Holtz man (2006) find that most 
of their galaxies have positive or zero gradient in [Mg/Fe] within 
the bulge-dominated region; as mentioned in the Introduction, they 
divide their sample in red (B~K > 4) and blue bulges, and, ac- 
cording to their analysis, the red bulges are overabundant in the 
centre, while the blue bulges have solar abundance ratios in the cen- 
tre, and the disc-dominated regions also have approximately solar 
abundance ratios. Therefore, our results are not necessarily in con- 
tradiction with theirs, as we will also see in Appendix IaI where we 
briefly discuss the radial profiles of the indices and the population 
parameters. 



5.2 Index - Hubble type relations 

The three panels in Fig.|9]present the relation between the central 
line indices (respectively, H/3, Fe5015 and Mgfo) and the morpho- 
logical type (see Table [Hi. Again, we use black symbols for the 18 
late-type spirals under investigation, red ones for the ellipticals and 
lenticulars and blue for the early-type spirals of the SAURON sur- 
vey. We see an overall increase of H/3 and a decrease of Fe5015 
and Mg& going towards later types. We also see that for ellipticals 
the range spanned in equivalent width is in general quite small, it 
becomes larger for lenticulars and Sa galaxies and smaller again in 
the later-type spirals. This result is illustrated quantitatively in Ta- 
ble [3] where we list for each index the range in equivalent width 
spanned in the three samples, together with the minimum and max- 
imum values assumed by the indices: as one can see, the range in 
H/3 is larger for the E/SO than for the Sa sample, while for the other 



indices the range is larger for the early-type spirals; in all three in- 
dices, the Sb-Sd galaxies span a narrower range than the Sa. In the 
same table, we also list the average value of the indices in the three 
samples, together with their standard deviation and the standard de- 
viation after subtraction of the observational errors, which indicates 
the amount of scatter that cannot be explained by the errors. 



5.3 Index - velocity dispersion relation 

Figure [To] shows the central Mg/? and H/3 indices against the cen- 
tral velocity dispersion of the stars (see Table]!}, in the left and right 
panels respectively. As in the previous Figures, the black symbols 
represent the late-type galaxies, the red the ellipticals and lenticu- 
lars and the blue the Sa observed with SAURON (central apertures). 
In these plots we express the indices in magnitudes, as it is often 
done in the literature, using the definition: 



index' — —2.5 x log ^1 — 



index \ 



(2) 



where index and index are measured respectively in A and 
magnitud es and AA is the width of the index bandpass (cfr. 
IPaper vj for the full index definition). 

A tight and well known relation between magne- 
sium indices (MgZ? or Mg2) and a hol d s for ellipti- 
cals ( see for example iTerlevich et alj ll98lL iBurstein et al.l 



| l988l iBender, Burstein & FabeJ Il99'3l|j0rgensen et alj Il99d 
Icolless et al." '1999*, Falcon-Barroso, Peletier & Balcell s 2002 
Jernardi et al. 2003, Worthey & Collobert 2003, Denicol o et al] 
I2OO5I) . In Fig. [To] (left panel) we overplot w ith a black solid line 
the relation obtained bv |j0rgensen et al. 1 ( Il996) for a sample of 207 
E and SO galaxies in 10 clusters: 



Mg2 = 0.196 X log(a) - 0.155, 



(3) 



after having converted the Mg2 index to Mgfe u sing a least-squares 
fit to afl the models from lVazdekis et alj ( Il996l) with Mg2 > 0.10, 
holding the relation: 



Mgfc = 0.489 X Mg2 + 0.014. 



(4) 



One could think that this transformation, using models that do not 
take into account differences in abundance ratios, could be mis- 
leading. But in practice the Mgi - Mg2 relation has very little 
scatter, even for different [Mg/Fe]. In any case, we also compare 
the position of our spirals in the (lo g (cr), Mgfe ) pla ne with the 
relations determined by Sanchez-Blazquez et al.l ( ]2006) in a fully 
empirical way for a sample of 98 early-type galaxies drawn from 
different environments (the field, poor groups, Virgo, Coma and 
some Abell clusters). The black dotted and dash-dotted lines over- 
plotted in Fig. [To] represent their fitted Mg& - o relation for low- 
and high-density e nvironments, respective ly, and lie both very close 
to the relation of lj0rgensen et al.l ( Il996l) . especially for high ve- 
locity dispersions; similarly, in the right panel we overplot their 
H/3 - a relations. The relations of J0rgensen et al. ( 1996) and 
ISanchez-Blazguez et alj ( I2OO6I) are both determined performing 
a fit down to vel ocity dispersions below 100 k m s~^ (~ 40 km 
s~^ in the case of lSanchez-Blazguez et alj2006l) , therefore explor- 
ing the range of velocity dispersions that we test with our late- 
type spirals. The small differences between these relations might 
be due to the different statistics and to the different environments: 
lj0rgensen et al.l i 19961) themselves noticed that the relation differs 
slightly from cluster to cluster, in the sense that galaxies in clusters 
with lower o have systematically lower Mg2 . We also performed a 
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Figure 9. Line-strength indices (in A) as function of morphological type. Black symbols are our sample galaxies (central aperture), red are the E and SO 
(central 1'.'335 aperture) and blue the Sa galaxies (central aperture) of the SAURON survey. The left panel refers to H/3, the middle one to Fe5015 and the right 
one to Mgh. Representative error bars are added at the bottom right of each panel; the black one refers to both the early- and late-type spirals samples, while 
the red one refers to the E/SO galaxies. 
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Sa 
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Fe5015range 
Mgbrange 


3.78 [1.13-4.90] 
2.50 [3.66-6.15] 
3.75 [1.51-5.26] 


3.10 [1.27-4.37] 
5.65 [0.31-5.96] 
4.05 [0.89-4.94] 


1.98 [1.96-3.94] 
3.71 [2.41-6.12] 
3.19 [1.20-4.39] 


(H/3) ± a 
(Fe5015) ± cr 
{Mgb) ± o- 


1.79 ±0.75 (0.67) 
5.45 ± 0.56 (0.41) 
4.33 ± 0.87 (0.77) 


2.39 ± 0.77 (0.57) 
4.51 ± 1.40(1.10) 
2.99 ± 1.07 (0.87) 


3.06 ± 0.56 (0.36) 
4.06 ± 1.04(0.74) 
2.19 ±0.77 (0.57) 



Table 3. Upper part: ranges spanned in equivalent width (A) by the three indices, for the three samples represented in Fig. |9l between brackets are indicated 
the minimum and maximum values. Lower part: average values of the indices in the three samples, with the relative standard deviations; the numbers in 
parenthesis indicate the standard deviation after subtraction of the observational error 



linear fit to our own data: the green solid line overplotted in both 
panels represents the corresponding relation. 

In Table |4] we list the coefficients defining the relations taken 
from the literature, together with the coefficients determined for our 
own sample, in the form: 

index' = aindex' + blndex' X log (cr). (5) 

As already pointed out by IPaper xj for the Sa galaxies of the 
SAURON survey (blue symbols in the Figure), it is clear from Fig. 
[TO]and from Table|4]that the Mgb - a relation for ellipticals repre- 
sents an upper envelope for the spirals; some of the spirals follow 
the relation (particularly if we consider the Sa), but the majority of 
them lie below it; in particular, the linear fi t to our own data (gree n 
line) runs almost parallel to the line from'j0r gensen et"al] j 19961) . 
but significantly below it. Similarly, in the case o f the H0 - a 
relation, the ellipticals represent a lower envelope jPaper x3) and 
many spirals lie above the relation. Deviations from the relations 
for early-type galaxies may be drive n by stochastic processes like 
star formation: as argumented by Schweizer & Seitzer (1992), the 
line defined by the galaxies of J0rgensen et al, (1996.) corresponds 
to old stellar populations, while the deviations would be due to 
younger stars. For our objects (black dots) the effect of young 
populations as deviation from the line defined by E/SOs galaxies is 
evident. 
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Table 4. Coefficients of the index - cr relation, in the form: index 
= a,; prfpT./ -I- bj„rff; T .' X log(cr); the first column lists the references 
(Jor= J0rgensen et aU 119961 : San^j=;Sanchez-Blazquez et al. 2006, low- 
density environment; San^fi3=:Sanchez-Blazquez et al. 2006, high-density 
environment; KG= this work), the second and third the coefficients for the 
relation with Mgfe , the fourth and fifth the coefficients for the relation with 
H/3'. 



6 AGE, METALLICITY, STAR FORIVIATION HISTORY 

The ultimate aim of measuring line-strengths is to understand the 
distribution of stellar ages, metallicities and abundance ratios, and 
to get hints on the star formation histories. Using a combination 
of line-strength indices, one can constrain stellar population mod- 
els, which allow to translate the line-strengths values into measure- 
ments of age, metallicity and abundance ratio. 

Trying to address this issue, we investigated two main ap- 
proaches: from one side, one can assume that a galaxy can be 
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Figure 10. Left panel: central Mgh , expressed in magnitudes, against central velocity dispersion, in units of decimal logarithm. The black symbols represent 
our own sample (central aperture), the red the E and SO and the blue the Sa galaxies of the SAURON survey. Right panel: central H/3 , expressed in magnitudes, 
against central velocity dispersion, in units of decimal logarithm; the colour code for the symbols is the same as in the previous case. The black solid line 
oveiplotted in the left panel represents the relation found by J0r2ensen et al. ( 1996) for early-type galaxies in 10 cluster The dotted and dash-dotted black lines 
overplotted in both panels are the relations obtained bv iSanchez-Blazquez et aU 12006) for low- and high-density environments, respectively, and the green 
solid lines are the relations determined using our own sample; see text in Section l531 for details. Representative error bars are added at the bottom right of each 
panel; the black one refers to both the early- and late-type spirals samples, while the red one refers to the E/SO galaxies. 



viewed as a Single Stellar Populations (SSPs), and compare obser- 
vations with theoretical SSPs in order to determine the best-fitting 
population; this is a 'classical' approach to stellar population stud- 
ies. From the other side, one can view a galaxy's population as the 
time evolution of an initial SSP; making some assumptions on the 
initial metallicity and the metal enrichment history, by comparison 
of observations and models one can retrieve information on the star 
formation history. Nowadays, this is a topic of greatest interest: re- 
cently, several groups have been working on it and developed tools. 
Here we mention noticeable examples such as STECMAP and its 
extension STECKMAP0 and MOPED. STECMAP and STECKMAP 
have bee n developed by Ocvirk and collaborators (for a descrip- 
tion, see lOcvirk et al.]|2006j. lOcvirk et alj|2006bl) . and try to re- 
cover the star formation history via spectral fitting, without any a 
priori assumption on the star formation history; the only condition 
imposed is that the star formation history is a smooth function of 
time; MOPED, developed by Panter and collaborators, is a method 
that allows to analyse huge quantities of spectra using the entire 
spectral range, and assumes that a galaxy can be viewed as super- 
position of single bursts (for details see Panter, Heavens & Jimenez 
2003). 

In this Section we will follow these two roads (SSP and recov- 
ery of star formation history), on the basis of our measured Lick 
indices. 



6.1 Single stellar population analysis 



Following the examples of IPaper VnSan dlPaperXll we used the 
single-burst stellar population models of iThomas et al.l ( |2003|) to 



http;//as tro.u-strasbg.fr/Obs/GALAXIES/stecmap_eng. html 



compare with our observations and estimate the ages, metallicities 
and abundance ratios of our galaxies. To avoid strong discretisation 
effects on the derived parameters, we interpolated the original grid 
of model indices, obtaining a cube of ~ 225000 individual models, 
with -2.25 s; Z < -1-0.67, 0. 1 ^ age 15 Gyr and -0.2 < [Mg/Fe] 
^ 0.5. We then determined the model closest to our observations 
(the set of three indices) for each Voronoi bin in our galaxies, via 
a minimisation technique, consisting in minimising contempo- 
raneously the distance between all of the observed and the model 
indices, weighted with the observational error on the line-strengths 
(see Section [3T| for the actual figures). We attributed to the galaxy 
bin the age, metallicity and abundance ratio of the selected model. 
The third row in Figures [6al6il shows the age, metallicity and abun- 
dance ratio maps obtained in this way. We are aware that it is an 
over-simplification to represent a galaxy's population using a SSP, 
however, we want to apply this method as a 'zero-order' estimate, 
before moving to somewhat more sophisticated ones. Fig.ll llplots 
against each other the age, metallicity and abundance ratio values 
calculated over the central aperture spectrum. Our sample galaxies 
are represented with black symbols. Overplotted with red symbols 
are the ellipticals and lenticular galaxies of the SAURON survey 
(Kuntschner et al., in preparation), while the blue symbols r epre- 
sent the early-type spirals of the SAURON survey jPaper xj) . De- 
spite the problems intrinsic to this approach, from this plot we can 
see that, as expected, our galaxies are on average younger and more 
metal poor than early-type galaxies, and have on average a value of 
[Mg/Fe] closer to zero, ranging from slightly sub-solar to slightly 
over-abundant. This confirms what we saw in Fig. [8] where we no- 
ticed that the late-type spirals were compatible with models with 
[Mg/Fe] =0. 

In order to assess the dependence of our results on the set of 
models used, we performed the following exercise. We took the 
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Figure 11. Central aperture values of the stellai' population parameters from our one-SSP analysis. Left panel: central age (in Gyr) against central metallicity 
(in decimal logarithm, with the solar metallicity as zero point); middle panel: central abundance ratio against central metallicity; right panel: central abundance 
ratio against central age. In black we represent our sample galaxies, while in red the E and SO and in blue the Sa galaxies observed with SAURON. 




Figure 12. Compari son between the ages (left panel, in Gyr) and metallicities (rig ht panel) obtain ed from our one-SSP approach using the models of 
Ixhomas et al.ll20O3l (horizontal axis, limitately to nearly-solai' abundance ratios) and IVazdekis|[T999l (vertical axis). The solid lines overplotted represent 
the 1 : 1 relation. 



model spectra from the SSP library of IVazdekij i ll 9991) and cal- 
culated the line-strength indices on them, obtaining a grid of 322 
model indices. Then we determined, galaxy by galaxy, the model 
that best matched our observations, deriving age and metallicity es- 
timates. Abunda nce ratios are fixed to solar values in the library of 
IVazdekislfl999l) . In Fig.[T2]we plot against each other the ages and 
metallicities obtained for the central I'.'S apertures of our galax- 
ies using the Thomas (horizontal axis) and the Vazdekis (vertical 
axis) models. Overplotted with solid lines are the 1:1 relations. We 
see that there are no dramatic differences between the two sets of 
measurements, though the Vazdekis models tend in several cases 
to predict slightly older ages and lower metallicities with respect 
to the Thomas ones. The small residual differences do not depend 
on the fixed abundance ratios of the Vazdekis models: the ages and 
metallicities on the horizontal axis in Fig. [12] are obtained for this 
test considering only the Thomas models closest to solar abundance 
ratios. The differences in the metallicity could be instead at least 
partially due to the coarseness in Z of the model grid obtained from 



the Vazdekis models. 

As also noticed in IPaper Xj. in some cases an unconstrained 
SSP fit cannot provide a suitable representation of a galaxy. In the 
case of a complex star formation history, not approximable with 
a single, instantaneous burst, the SSP-equivalent parameters must 
be interpreted with caution, as a 'zero-order' estimate, as we al- 
ready stressed. As an example, we remind the reader the case of th e 
star formation ring in NGC4321 discussed byj Allard et al. | |2006|). 
They show that in the ring the line-strength indices suggest the su- 
perposition of two components: a young stellar population and an 
old and metal-rich one. When they force a si ngle SSP, they obtain 
low metallicity and old age. But it is known iZaritskv at alj|l994h 
that the HII regions of this ga laxy have super-sol ar metallicity: the 
one-SSP solution obtained bv lAllard et alj ( l2006h represents there- 
fore an inconsistent description. As clearly visible from the first 
panel in Fig. [TT] some of our galaxies also populate the region of 
the (metallicity, age) plane characterised by low metallicities (be- 
low — 1.0) and old ages (above 10 Gyr). This is the case for the 
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central aperture of NGC 864, for which our technique selects a 
model with Z ~ —1.14 and age ~ 11 Gyr. The line-strength maps 
for this object (see Fig.l6bt show that in the very center the Fe5015 
and Mgfc indices assume low values, and so does H/3, on a large 
inner region. The low indices drive the metal poor old population, 
but a SSP description might be too unrealistic for this object. 

More generally, if for example a galaxy has undergone two 
separate bursts of star formation, the SSP-equivalent age will be bi- 
ased towards the age of the youngest stars, and the SSP-equivalent 
metallicity will be biased towards the met allicity of the old popu la- 
tion, as studied and described in detail bv lSerra & Trageil ( l2007t) . 

In order to put some physically motivated constraints to our 
one-SSP analysis, we decided to constrain the metallicity of our 
objects in a narrow range around the value in the model grid clos- 
est to the one determined from a relation between central stellar 
velocity dispersion and metallicity for early-type galaxies. In prac- 
tice, we aim to apply a scaling relation obtained for spheroids (el- 
lipticals and lenticulars). Our late-type galaxies do not host rele- 
vant bulge components (Ganda et al., in preparation), so the re- 
gions that we probe with our data do not fall under this category. 
Also, unlike for spheroids, o is not a good mass indicator in disc 
galaxies. But in any case, using a 'cr relation' for spheroids is a 
'zero-order' approach to extrapolate the metal content of the old 
stars in low-mass galaxies. We investigated several possibilities for 
a metallicity-sigma relation for early-type galaxies. The resulting 
relations can differ from each other significantly, particularly at the 
low-o" end, because of differences in the methods, in the models 
used, and in the underlying assumptions. 

A possible and common way is to go through stellar popula- 
tion models, estimate the ages and metallicities by comparison of 
observations and SSP models, and then fit a relation with o. This 
approach has been extensively investigated in the literature; we will 
refer to it as ' the model appr oach' . As a prototype, we refer to the 
work of iThomas et~^ feOOSi) . who studied 124 early-type galax- 
ies in high- and low-density environments and derived their ages, 
metallicities and element ratios. For the galaxies in high-density 
environments, they retrieved the following relation: 

Z = 0.55 X log(o-) - 1.06. (6) 

See also lKuntschned fcOOOl) for a similar result. 

Another possibility, to which we will refer as 'the index - a 
approach', is to assume a set of index - a relations for early- type 
galaxies, extrapolate them to the range of a under investigation, 
from the measured a infer the indices and from these, via a com- 
parison with models, an estimate for the metallicity. In practice, we 
used the H/3 - cr, Fe5015 - cr , Mgfe - a relations published in 
IPaper Vll for the 48 ellipticals and lenticulars of the main SAURON 
survey. Using these relations and the measured values for the cen- 
tral aperture velocity dispersions of our galaxies, we obtained a 
set of index values. We then compared them with the models of 
iThomas et al.l ( l2003h via the previously described minimisation, 
but choosing only among models with solar abundance ratios and 
age « 12.6 Gyr. In this way we come to the estimate of the metal- 
licity of an old spheroid at the observed velocity dispersion. We 
will adopt the estimates obtained in this way to put a constraint on 
the metallicity. 

Another approach, which does not involve any stellar popula- 
tion modeling, relies on empirical calibrations and tight observed 
scaling relations between metal indices and velocity dispersion. 
Consider the Mg2 - cr relation fo r elliptical and lenticular galax- 
ies given bv ij0rgensen et al.]| 19961 (Equation [Sj. If we combine this 
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Figure 13. Relation between central metallicity and central velocity disper- 
sion (in decimal logarithm and units of km s~^) for the E and S O galax- 
ies (Kuntschner et al., in preparation, red symbols) and the Sa I Paper Xl 
blue symbols) of the SAURON survey and for our late-type spirals (black 
dots), from our one-SSP analysis. The overplotted lines represent the rela- 
tion metallicity - velocity dispersion obtained in different ways, drawn on 
the cr range of our data only: the dotted line reproduces Equation|6]('model 
approach'), the dashed line Equation[8]('empirical approach') and the solid 
lines represents a relation based on the assumption of a set of index - cr 
relations. See text for a more complete description. 

relation with the Mg2 - metallicity calibration from Buzzoni, Gari- 
boldi & Mantegazza (I992f|: 

Mg2 = 0.135 X ^ + 0.28, (7) 

we obtain a relation between metallicity and velocity dispersion, 
holding for early-type galaxies l3[ 

Z = 1.452 X log(cr) - 3.222. (8) 

We call this last approach 'the empirical approach'. 

In Fig. [13] we show in the (log (cr), Z) plane the distribu- 
tion of the elliptical and lenticular galaxies in the SAURON sur- 
vey (red symbols, from Kuntschner et al., in prepa ration), t he Sa 
galaxies from the same survey (blue symbols, from lPaperXlh , to- 
gether with our own galaxies (black dots). All the values refer to 
the central apertures, and the metallicities are obtained from the 
one-SSP analysis. The lines overplotted in the Figure represent the 
various metallicity-velocity dispersion relations for E/SO galaxies 
described above. The dotted line reproduces Equation [6] (from 'the 
model approach'), the dashed line Equation [8] (from 'the empirical 
approach') and the solid line connects the points representing our 
galaxies in the (log (cr), Z) plane according to 'the index - cr' ap- 
proach, the second one we presented. 

We re-evaluated our age estimates by imposing a constraint 
on the metallicity, allowing it to vary only in a narrow range (± 0. 1 
dex) around the metallicity given, galaxy by galaxy, by the solid 

^ These authors calibrate the dependence of Mg2 from stellar parameters 
using an observed stellar library and find a general relation, that reduces 
to the one reported here (Equation [7) in the parameter range spanned by 
elliptical galaxies. 

^ We warn the reader about the fact that the Mg2 index has a slightly dif- 
ferent definition in the papers of |j0rgensen et alj )l99d) and lBuzzoni et alj 
1 1992), and that we have not taken this difference into account in deriving 
Equation[8] 
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Figure 14. Same plots as in Figure [T7I but referring to the case with constrained metallicity (see text in Section ISTI for more details), for our late-type spirals 
(black symbols); for the E/SO (red symbols) and the Sa (blue symbols) galaxies no constraint on metallicity is applied. 




age free (GvO 

Figure 15. Constrained (y-axis) vs unconstrained (x-axis) age (in Gyr), 
measured on the central apertures of our galaxies with the one-SSP ap- 
proach. 



black line in Fig. [T3] ('the index - o approach' discussed above). 
We preferred this approach to the other two because it is based 
on empirical relations, and does not rely exclusively on models, 
unlike 'the model approach', and because it uses more than just 
magnesium to establish the metallicity of the old population, con- 
trary to 'the empirical approach'. Fig. [14] is equivalent to Fig. II II 
but with the constraints on metallicity. In Fig. [15] we plot the con- 
strained ages against the unconstrained ones: we can notice that for 
very young ages the values are close, and differ more and more go- 
ing to older ages. The fact that in some cases the constrained age 
differs from the unconstrained might indicate either that the im- 
posed metallicity is inappropriate or that the one-SSP assumption is 
wrong, and the galaxy's star formation history is not approximable 
with an instantaneous burst. 
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Table 5. Population parameters from our one-SSP analysis. The subscript 
'c' refers to the quantities obtained in the case where we constrain the metal- 
licity within a narrow range (± 0.1 dex) ai'ound the metallicity given by the 
Z - (T relation. Age is in Gyr and metallicity in decimal logarithm, with solar 
metallicity as zero point. 



In Table|5]we list the population parameters measured, both in 
the unconstrained and in the constrained case, on the central aper- 
ture of our galaxies. 



6.2 Two-population analysis 

A galaxy is most likely characterised by stellar populations that are 
more complex than SSPs; therefore, in order to try to build a more 
realistic picture of our galaxies, we investigated more sophisticated 
approaches. As an attempt, we explored a two-SSP scenario, de- 
scribing the galaxy as the superposition of an old population and a 
younger one. In practice, we want to describe the galaxy population 
as: 

popjtot = (1 — /) X pop.old + / X popjyoung , (9) 
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Figure 16. Fitted parameters from our two-SSP analysis, measured on the central apertures. In the left panel we have the age (in Gyr) against the metallicity 
(in decimal logarithm, with the solai' metallicity as zero point) of the young component popjyoung; in the middle its mass fraction against its metallicity; in 
the right panel, the mass fraction of the young component against its age. 



where pop_old is a ^ 12.6 Gyr old SSP with metalhcity fixed to 
the value close to the one obtained from the Z - cr relation (on the 
basis of what we called 'the index - a approach') and popjyoung 
is a SSP with age below 5 Gyr and free metallicity; the limit of 5 
Gyr was chosen in order to have a population younger than the 
old component, but at the same time not limited to very young 
ages; / is the mass-fraction of the young component and varies be- 
tween and 1 . To do this, we take the model spectra from Vazdekis 
( Il999f) and select, galaxy by galaxy, the 12.6 Gyr old model with 
the metallicity closest to the one given by the solid line in Fig. 1131 
among those available in the model grid, and combine it with all 
the other SSP models in the library, with varying mass fractions. In 
this way we obtain for each galaxy a set of 6762 'composite model 
spectra', with varying age and metallicity of the young component 
and relative mass fraction. On these spectra we compute the line- 
strength indices, obtaining for each galaxy a g rid of 'composit e 
model indices'. Since we use the model library of lVazdekisI jl999l) . 
the abundance ratio is fixed to solar values, which for our objects 
is a quite fair approximation, as we saw in the previous Section. 
We then use the same technique as for the one-SSP analysis and 
determine, spectrum by spectrum, the composite model closest to 
our observed line-strengths. This gives us an estimate for the age 
and metallicity of the young component and its mass fraction. 

Figure [T6]plots against each other the parameters fitted on the 
spectrum averaged in a central aperture of l'.'5 radius. We see that 
the young component spans a range in age from below 1 to ~ 4.5 
Gyr, a range in Z from — 1 .38 to -1-0.2 and that its mass contribution 
to the global galaxy population / varies from ~ 10% (NGC2805) 
to ^ 95% (NGC 1042, 2964, 3949, 4487). 

The main problem of this approach is the degeneracy of the 
parameter space, which must act as a warning for the reader against 
the robustness of the ages and mass fraction estimates for the young 
component presented in Fig. [16] This degeneracy is illustrated in 
Figures [it] and [TS] In Fig. [17] we plot, for different choices of the 
metallicity of the old component and solar metallicity for the young 
one, the points representing the models in the ([MgFe50], H/3) di- 
agram. Different colours represent different weights of the young 
component in the composite models, from O.I (black) to 0.9 (gold). 
From these plots it is evident that it is not possible to significantly 



discriminate the effects of age and mass fraction of the young com- 
ponent. Fig.[T8]shows the projection of the Ax^ space (metallicity, 
age and mass fraction of the young component) on the (age, frac- 
tion) plane and the contours representing the 1,2,3-cr confidence 
level, based on the Ax^ for two degree of freedom, for the central 
aperture of some of the galaxies. We see that very large regions of 
the parameter space are equivalent in terms of and that even 
in cases where the age is rather well constrained (see NGC 2805), 
the mass fraction is badly determined. Therefore, we see that we 
cannot extract solid conclusions from this method, despite its inter- 
esting principle, and we decided not to push it further and not to 
present the two-dimensional maps of the inferred parameters. 



6.3 Continuous star formation 

'The star formation history of galaxies is imprinted in their inte- 
grated light', state Bruzual & Chariot in a paper where they present 
a new model for the spectral evolution of stellar populations at 
ages between 1 x 10^ and 2 x 10^*^ yr for a range of metallici- 
ties teruzual & Charlolll2003h . In practice, they use the isochrone 
synthesis technique to compute the spectral evolution of stellar pop- 
ulations. This technique relies on the fact that any stellar population 
can be expanded in a series of instantaneous starbursts or SSPs. So, 
the spectral energy distributi on at time t of a stellar pop ulation can 
be written as in Equation I in iBruzual & CharlotI bOOSi) : 

Fx(t)^ I ^{t-t')y.Sx[t',(,{t-t')\dt' , (10) 

where ^(f) and C,{t) are respectively the star formation rate and the 
metal-enrichment law, as a function of time, and Sx \t' , C,{t — t')] 
is the power radiated by an SSP of age t' and metallicity ({ t — t') 
per unit wavelength and unit initial mass. Bruzual & Chariot 
compute their models with various choices for the initial mass func- 
tion IMF, the metallicity of the initial SSP and the stellar evolu- 
tion prescription. They also provide as online material some tools 
that allow the user to compute the spectral evolution -in time- of a 
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Figure 17. Composite models for different choices of the metallicity of the old component (indicated at the top of each panel) and solar metallicity of the 
young one represented by points in the ([MgFe50], H/3) plane. Different colours refer to different values of the mass fraction of the young component in the 
composite model. The effects of age and fraction of the young component on the line-strength indices of the composite models are degenerate. 
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Figure 18. Contours of Ax^ projected on the plane (age, fraction) of the young component at the 1,2,3-(t confidence level, for the central apertures of 
NGC772, NGC2805, NGC4254; the red symbol marks the best-fitting model. These plots cleai'ly show that the two-SSP approach is unable to constrain the 
fitted parameters, and that age and mass fraction of the young component are degenerate. 



SSP for different star formation history scenario^ We chose to use 
these tools in order to apply an evo lutionary analysis to our data. 
Starting from the SSP models from lBmzual & CharlotI Ew^ ob- 
tained using the Chabrier (2003) IMF, for different choices of the 
initial metalhcity (namely, Z = -1.69, -0.69, -0.39, 0.0, 0.39), 
we computed the time evolution with constant and with exponen- 
tially declining star formation rate: respectively, ip{t) — const 
and ip{t) = [IMq + eMpG{t)]T~'^exp{-t/T), where r is the 
e-folding time-scale and Ad pc (i) is the mass of gas processed into 
stars and re-ejected into the interstellar medium at the time t; the 
parameter e sets the fraction of the ejected gas that is re-used in 
new star formation episodes; we only explored the cases e = (no 
recycling) and e — 1 (all of the ejecta go into new stars). In any 
case, the actual value of e does not change our results significantly. 
These are both plausible scenarios: the constant star formation his- 



" The program we used are available from the web page 
[http://www.cida. ve/'^bruzual/bc2003| . 



tory has long been considered as a suitable scenario for spirals: 
'star formation has proceeded at relativel y constant rate o ver the 
lifetimes of most late-type spiral galaxies' ( lKermicut3l I983I) . while 
the exponentially declining star formation describes the situation of 
a starburst of duration r; we notice that the case r — > cxa approx- 
imates the constant star formation scenario, while the case r 
reproduces an instantaneous burst or SSP. We computed the spec- 
tral evolution for the chosen metallicities, for both cases e = and 
e = 1 and for several choices of r: r = 0.5, 1.0, 3.0, 5.0, 7.0, 10.0, 
12.0 and 15.0 Gyr. The tools provided by Bruzual & Chariot also 
include an option for computing the effect of attenuation by dust, 
but we decided not to use this possibility, in order not to further 
increase the number of parameters at play. 

For all of the selected cases we calculated the time evolution 
of the Lick indices using the programs from Bruzual & Chariot and 
focussed on the situation at t = 10 Gyr, building a model grid of 
Lick indices for different metallicities, r and e (0 or 1) and age fixed 
at 10 Gyr and then compared them with the line-strength indices of 
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Figure 19. Central aperture values for the e-folding time-scale r against central velocity dispersion cr (both in units of decimal logarithm), in an exponentially 
declining star formation scenario; the t values are obtained selecting among models with age = 10 Gyr and e = 1; the error bars come from a rough estimate 
on the basis of the l-cr contour plots of the projection of Ax^ on the (r-metallicity) plane. 



the SAURON spectra, via the minimisation technique used also 
in the previous Sections. In this way, we come to an estimate for 
the time-scale r, e and the metallicity. Fixing the age to 10 Gyr is 
equivalent to assume that the galaxy existed for 10 Gyr, and to try 
to describe how stars were formed over that time. We tested that the 
fitted parameters (r and metallicity) had negligible variations when 
fixing the value of e, therefore we repeated our estimates allowing 
only models with e = 1 and will refer to this case from now on. 

In Figures [6a|6i| the top right panel shows, galaxy by galaxy, 
the two-dimensional maps of r, obtained for the exponentially de- 
clining SFR, for age fixed to 10 Gyr and e = 1. Interesting features 
of the individual r maps, whenever present, are highlighted in Ap- 
pendix |B] As for the metallicity maps, obtained in the minimi- 
sation and not shown in this Paper, we can state that in the majority 
of cases they resemble in their structure the metallicity maps ob- 
tained in the one-SSP approach; the actual numbers are different, 
since the models used in the two approaches are different and also 
the metallicity ranges and values explored by the model grids dif- 
fer, but the spatial structure is generally similar. 

If we focus on the values obtained for the central apertures, 
we find a correlation between the best-fitting r and the central ve- 
locity dispersion, as shown in Fig. [19] where we plot log(T) against 
log(o-). We see that higher a galaxies, which are larger and more 
massive, tend to have shorter star formation time-scales, more con- 
sistent with an instantaneous burst scenario, while at the low-a end 
the smaller galaxies tend to have larger r, indicating a star forma- 
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Table 6. Best-fitting values for the star formation time-scale r (in Gyr) 
measured on the central apertures of our galaxies, obtained for the choices 
age =10 Gyr and e =1, in an exponentially declining star formation rate 
scenario. 

tion history spread over time. Therefore, the time-scale r represents 
a physically meaningful parameter, parametrizing a smooth transi- 
tion from SSP-like to constant SFR. 

Figure |20] shows that rather tight correlations exist also be- 
tween the central r and the central line-strength indices. These 
plots tell us that old, metal-rich galaxies (low H/9, high Fe5015 and 
Mgb values) have shorter star formation time-scales, while young 
and more metal poor galaxies (higher H/3, lower Fe5015 and Mgb) 
have a star formation history more extended in time: we find again 
a smooth transition between a SSP-like and constant star formation. 
The central aperture values of r used in Figures[T9]and|20]are listed 
in Table |6l 
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Figure 20. Central aperture values for the e-folding time-scale t (in Gyr) of an exponentially declining star formation history, against the central line-strength 
indices (expressed in A): H/3 in the left panel, Fe5015 in the middle and Mgb in the right one; the r values are obtained selecting, in the fit, among models 
with age = 10 Gyr and e = 1 (see text for details). 



We also investigated the constant star formation history sce- 
nario, by selecting the corresponding option when using the pro- 
grams from Bmzual & Chariot, building a grid with the model in- 
dices for the exponentially declining star formation rate and adding 
those for constant star formation rate, for the various metallicity 
values selected and for e = 1 and age fixed to 10 Gyr. When ap- 
plying our X minimisation procedure, we found that the constant 
star formation was preferred in the regions where the fit with the 
exponential star formation only had selected the highest values al- 
lowed for the time-scale r. This confirms our interpretation of r as 
a parameter describing a continuous transition from constant SFR 
to SSP. 

We then repeated the analysis for different values of the fixed 
age of the evolved models, spanning the range 1-18 Gyr and found 
that the r values depend on the fixed age, given the existence of a 
strong degeneracy between age and r, while the metallicities are 
more robust. 

The main interpretative difficulty of this method is indeed due 
to the degeneracies in the parameter space, which are difficult to 
break with the few observed Lick indices that we measure. Nev- 
ertheless, by inspecting the contours of Ax^ in the (r-metallicity) 
plane, for several choices of the fixed age of the evolved models, we 
can in some cases exclude certain star formation scenarios in favour 
of others, even though we cannot give a precise, quantitative esti- 
mate of the fitted parameters. For example, in the case of NGC 488, 
if we look at the two-dimensional Ax^ space (metallicity, r), for 
a fixed age of the evolved models), and look at the contours in the 
(r-metallicity) plane representing the 3-a confidence level, based 
on the Ax^ for one degree of freedom, we notice that it lies in a 
region far from the high-r end, and this happens independently on 
the chosen values of the fixed age. Even if we perform the min- 
imisation by allowing free age, r and metallicity, and examine the 
projection of the Ax^ contours in the (r-age) plane, we see that for 
all ages we can exclude large r values at a high confidence level. 
This must be regarded as an exercise, since we are not interested 
in fitting simultaneously age and r, but it gives us insight on the 
impact of the degeneracies and of our assumptions on the robust- 
ness of our conclusions. Thus, for NGC 488 we can safely state 
that we can exclude a constant star formation scenario in favour of 
an exponentially declining star formation rate, with a short time- 



scale T (closer to an instantaneous burst than to a constant SFR). 
This is illustrated in the first column in Fig. I21al Similar hints for 
the star formation history of other galaxies will be given, whenever 
possible, in Appendix |B] dedicated to the description of the single 
galaxies individually. For all galaxies, the Ax^ contours for the fit 
at different fixed ages and with free age are presented in Figures 
|21a|21c| 



7 SUMMARY 

In this Paper we presented the first integral-field measurements of 
line-strength indices and stellar population parameters for a sam- 
ple of late-type spiral galaxies. This work is based on data of 18 
spirals spanning a range in type between Sb and Sd, observed with 
SAURON, which provides measurements of the Lick indices H/3, 
Fe5015 and Mgb over a field of view covering the nuclear region 
of our objects. 

Here we summarise the main achievements of the Paper: 

• We presented the two-dimensional maps of the line-strength 
indices and analysed the values extracted on a central aperture of 
1'.'5 radius. 

• We looked at the relations with morphological type: we found 
that late-type spirals do differ from earlier-type galaxies in terms 
of their populations, having globally higher H/3 values and lower 
metal indices than ellipticals, lenticulars and early-type spirals. 
This suggests younger ages and lower metallicities. 

• We investigated the relations between indices and central ve- 
locity dispersion: late-type spirals do not obey the index - a rela- 
tions found for early-type galaxies: our galaxies lie above the H/3 
- a and below the Mgb - a relations determined for the elliptical 
galaxies, that correspond to old populations. 

• We determined SSP-equivalent ages, metallicities and abun- 
dance ratios by comparison of our observed indices with models: 
ages are mainly young, metallicities low and abundance ratios close 
to solar; the range spanned in the parameters is narrower than for 
early-type spirals. 

• We investigated the possibility of recovering the star forma- 
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Figure 21a. First column: projections of the 1-, 2- and 3-(t Ax^ contours on the (r-metaUicity) plane for the central aperture of NGC 488, for different fixed 
ages in the fit (5 Gyr, first row, 10 Gyr, second row, 18 Gyr, third row); projections of the 1-, 2- and 3-(T Ax^ contours on the (r-age) plane for the central 
aperture of NGC 488, for the fit performed with free age, r and metallicity, fourth row; second to fifth column: as in the first column, for NGC 628, 772, 864, 
1042 and 2805; in all plots, the red symbol indicates the best-fitting model. These plots demonstrate that for NGC 488, 628 and 772 we can exclude at a high 
confidence level a constant star formation scenario, while, for example, for NGC 2805 high values of t are not highly unlikely at all ages. 
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Figure 21b. As in Figure|2Ti tor NGC 2964, 3346, 3423, 3949, 4030 and 4102. 
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Figure 21c. As in Figure[2Ta] for NGC4254, 4487, 4775, 5585, 5668 and 5678. 
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tion history: we assumed that the galaxies formed 10 Gyr ago, 
and that they experienced an exponentially declining or constant 
star formation rate since then, and compared our indices with the 
evolved models of iBruzual & C hariot (2003); the fit returned an 
estimate for the e-folding time-scale for star formation t, which 
parametrizes a smooth transition from a SSP-like star formation 
history to a constant SFR: small r approximate an instantaneous 
burst and very large r reproduce a constant SFR; for the central 
apertures of our galaxies, we found a large range in r values, from 
1 to 15 Gyr, covering basically the whole parameter interval that 
we explored. 

• Interestingly, we found a trend between the fitted r and the 
central velocity dispersion, in the sense that galaxies with large a 
tend to have small r and galaxies with low a tend to have large r. 

Our observations and results nicely fit in the scenario summarised 
by Kennicutt ( 1998): there are two modes of star formation; one 
takes place in the extended discs of spiral and irregular galaxies, 
the other in compact gas discs in the nuclear regions of galaxies. 
The first has been detected mainly via Hq surveys and it is strongly 
dependent on the morphological type: the star formation rate in- 
creases by a factor of ~ 20 going from Sa to Sc galaxies. Similar 
information is derived from observations of the UV continuum and 
broadband visible colours. The second mode of star formation can 
be found in the circumnuclear regions; the physical conditions in 
these regions are in many cases distinct from the more extended star 
forming discs. The circumnuclear star formation is characterised 
by the absolute range in SFR, much h igher spatial con centration of 
gas and stars and its burst-like nature jKeiuiicutjl 1998 ^. Circumnu- 
clear star formation is largely decoupled from Hubble type. Studies 
on the dependence of nuclear Hq emission in star formi ng nuclei 
as a function of galaxy type ( IStauffeilll98"l iKeelll 19831 Ho, Fil- 
ippenko, & Sargent 1997) showed that the detection frequency of 
HIT region nuclei is a monotonic function of type, increasing from 
0% in elliptical galaxies to 8% in SO, 22% in Sa, 51% in Sb, and 
80% in Sc - Im galaxies. Among the galaxies with nuclear star for- 
mation, the Hq luminosities show the opposite trend; the average 
extinction-corrected luminosity of HII region nuclei in SO - Sbc 
galaxies is nine times higher tha n in Sc galaxies. T he conclusion of 
these studies, as summarised bv lKennicuttI ( Il998h . is that the bulk 
of the total nuclear star formation in galaxies is weighted towards 
the earlier Hubble types, even though the frequency of occurrence 
is hig her in the late types. According to iKormendv & Kennicutll 
( 12004 ) ■ the central star formation accounts for 10 - 100% of the 
total SFR of spiral galaxies. The highest fractions occur in early- 
type s piral galaxies, which typ ically have low SFRs in their outer 
discs iKenmcutt & Kentll 19831) . 

The work presented in this Paper supports this picture: if we 
put together the information gathered from the sample of E/SOs of 
Paper VI, the sample of Sa galaxies of Paper XI and our own sam- 
ple, and look at the behaviour along the Hubble sequence, we find 
that early-type spirals show a larger range in age than both E/SOs 
and late- type spirals; this is consistent with a scenario where the 
star formation in elliptical and lenticular galaxies is not very im- 
portant, in early-type spirals it is dominated by short bursts and in 
late-type galaxies it is more quiescent. This is also supported by 
our finding that big spirals tend to have an SSP-like star formation 
history, while for smaller ones the star formation history is better 
approximated by a constant over time. 

The following steps in the analysis of the SAURON data here 
presented and discussed will be a full investigation of the spatial 
information, with a more detailed study of the radial profiles (here 



shown in Appendix O and the gradients, and an interpretation of 
the connection between features in the line-strengths maps and in 
the morphological appearance. The interdependencies between the 
populations and the kinematical structures ( Ganda et a l. 2006) will 
also be examined. 

The main limitation of our work is the degeneracy in the 
parameter space in the different approaches we investigated. Our 
analysis relies on three data-points only (the line-strength indices 
within the SAURON spectral range), which are too few to break 
the degeneracies and constrain the fitted parameters. Therefore, it 
is urgent to complement the data presented here with spectra with 
a broader spectral coverage. Observations at other wavelengths - 
both shorter and longer- would help building a complete picture: 
both the UV and the NIR spectral regions can provide popula- 
tion indicators much less affected by the model degeneracies found 
in the optical. In addition, different sub-populations contribute to 
the integrated light in different spectral ranges, therefore a multi- 
wavelength study would be a step forward towards a full under- 
standing of the formation and evolution of our objects. 
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APPENDIX A: RADIAL BEHAVIOUR 

We briefly investigated the radial variations of the line-strength in- 
dices and age, metallicity and abundance ratios. From the corre- 
sponding maps, we calculated radial profiles of H/3, Fe5015, Mgfc 
and age, Z and [Mg/Fe] from the one-SSP approach, by averaging 
the maps on concentric annuli of (/.'S width. We present these pro- 
files in Figures lAlalAlcl These plots show that there is a wide va- 
riety of radial behaviours among late-type spirals: the line-strength 
indices can be roughly constant in radius (H/3 profile for NGC 488), 
rising outwards (even very steeply, see for example the H/9 profile 
for NGC 5678), rising inwards (many Fe5015 and Mg& profiles). 
In some cases there are also complicated structures, with 'oscil- 
lating' profiles: see for example the Fe5015 profile for NGC 864 
or NGC 2805 and several others. Some of the 'oscillating' profiles 
might just reflect the noisy nature of the data for some low-surface 
brightness objects (for example: the Fe5015 profile for NGC 2805), 
but others are probably related to real variations of the indices 
within the galaxy. Despite all this variety, we conclude that in 
most cases the H/3 absorption strength increases with radius, while 
Fe5015 and MgZ? decrease. The fact that the metal lines show most 
often negative gradients was already known for early-type galax- 
ies, while in the case of H/3 for elliptical galaxies the gradients are 
consistent with a flat relation with radiu s, or present a mild positive 
outwards radial gradient (see lPaper Vll and references therein). 

The stellar population parameters show the same heteroge- 
neous behaviours. Roughly, the age shows most often a complex 
behaviour with radius, reflecting features recognizable in the maps 
in Figures I6al6il whil e the metallicity in m any cases is decreas- 
ing, in agreement with lMoorthv & HoltzmanI l i2006l) , who measure 
mostly negative metallicity gradients for a sample of SO-Sc galax- 
ies; in the majority of the galaxies, the metallicity profile resem- 
bles closely the aspect of the Fe5015 profile. Looking at the abun- 
dance ratio profiles in more detail, we see that the central values 
are usually very close to solar; many of the profiles stay approxi- 
mately constant with radius, after an initial increase in the bulge- 
dominated region, which can be very tiny (for a quantitative de- 
scription of the bulges in these galaxies, we refer the reader to 
Ganda et al., in preparation). This i s also consistent with the find- 
ings of lMoorthv & Holtzma j ( |2006|) . that we briefly summarised in 
Section fSTI they suggest that the bulges have positive or null gradi- 
ents in abundance ratio with increasing radius, while the discs have 
solar abundance ratios. In early-type galaxies, instead, the radial 
profiles of abundance ratios are rather flat, as found by Kuntschner 
et al. (paper in preparation) for the SAURON E/SO galaxies and by 
ISanchez-Blazguez et al.l ( l2007h for a sample of 1 1 E/SO galaxies; 
the latter authors do find radial abundance ratio gradients, both pos- 
itive and negative, but rather shallow. 



APPENDIX B: INDIVIDUAL GALAXIES 
NGC 488 

NGC 488 is an unbarred Sb galaxy and it is the most elliptical-like 
(featureless) object within our sample. NGC 488 has regular 
stellar and gaseous kinematics, with regu lar rotation and ve locity 
dispersion smoothly decreasing outwards jGanda et al.l l2006!). The 
line-strength maps present high values for Fe5015 and Mgfo and 
low values for H/3 over the central regions. As one can notice from 
the values reported in Table |2l this galaxy has the highest Fe5015 
and Mgfc in the sample, and the lowest H/3. The ages are mainly 
old and the metallicities high. This is indeed the most metal-rich 



galaxy in the sample (the central metallicity from the one-SSP 
approach is ~ 0.49, then it decreases outwards). This galaxy is 
well described by a one-SSP: as shown in Fig. |21a| (first column) 
on the basis of our continuous star formation analysis, NGC 488 
is very unlikely to have experienced a constant star formation 
over time; the observed indices are instead reproduced by a short 
starburst, approximating a SSP. The time-scale of star formation r 
is short over the whole field. From the profile in Fig. lAlal we can 
see that the [Mg/Fe] profile is approximatively constant around the 
value ~ 0.2: abundance ratios are slightly super-solar, confirming 
the fact that this galaxy resembles an elliptical more than the later 
objects in our sample. Other suggestions in this sense come from 
the position it occupies in the diagrams in FigureslTllSllllI where 
it lies in the region occupied by the ellipticals. 



NGC 628 

NGC 628, known a lso as M74, is a grand-design Sc galaxy. 
ICornett et ai] ( Il994l) concluded that the star formation history of 
NGC 628 varies with galactocentric distance; Natali, Pedichini & 
Righini (1992) suggested that the galaxy could be seen as an inner 
and an outer disc characterised by different stellar populations; 
according to them, the transition between the two regions is located 
at ~ 8-10 kpc from the centre. Our observations of NGC 628 are 
disturbed by the presence of a foreground star which falls close to 
the centre (« 13" southern). The stellar kinematics is characterised 
by slow projected rotation, and velocity dispersion decreasing in 
the c entral zones, indic ating a cold central region, maybe an inner 
disc l lGanda et al.l2006l) . This galaxy is part of the 'control s ample ' 
of non-active objects studied by ICid-Femandes et alj 1 20041) . 
iGonzalez Delgado et all ( |2004|) . Icid-Femandes et al.l ( l2005h . who 
analyse the nuclear spectrum, describing it as a mixture of 10^-10^ 
yr-old and older stars. 

As can be seen from the profiles reported in Fig. lAlal 
NGC 628 has ages (one-SSP approach) below 3 Gyr out to « 15", 
and becomes older at larger radius; the features are roughly in 
agreement with the above mentioned findings of Cid-Femandes 
and collaborators; the same trend in age is obse rved, on the 
common radial range, also bv lMacArthur et al] ( l2007h . on the basis 
of CMOS data independently analysed. The metallicity decreases 
outwards. The 'bumps' recognizable in the Fe5015, Mgfe and (but 
less pronounced) metallicity profiles at ~ 13" from the centre are 
due to the foreground star. From our continuous star formation 
analysis, as shown in Fig. |21a| (second column), we can state that 
this galaxy is not described by a constant star formation scenario. 



NGC 772 

NGC 772, also called Arp78, is an Sb galaxy characterised by a 
particularly strong spiral arm; it forms a pair at 3f3 with the E3 
galaxy NGC 770. From the kinematical point of view, NGC 772 is 
characterised by a clear drop of the stellar velocity dispersion in the 
central zones, possibly corresponding to an inner disc, as suggested 
also by the a nticorrelation betw een velocity and the Gauss-Hermite 
moment h3 jGanda et ai] |2006b . In the line-strength maps we can 
clearly see an extended central region with high Fe5015 and Mg/? 
and low H/3; indeed, from the central aperture values reported in 
Table|2]we can see that this galaxy is one of those with the highest 
metal indices, surpassed only by NGC 488. In the H/3 map it is 
possible to see, surrounding the central depression, regions with 
higher values corresponding to the spiral arms. A similar behaviour 
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Figure Ala. Azimuthally averaged radial profiles of the line-strength indices H/3, Fe5015, Mgfo, (expressed as equivalent widths and measured in A) and 
of the one-SSP analysis parameters age (in Gyr), metallicity and abundance ratio against galactocentric distance (expressed in arcsec); every row presents a 
different galaxy. See text in Appendix |A]for a more detailed description. Here we show NGC 488, NGC 628, NGC 772, NGC 864, NGC 1042 and NGC 2805. 



can be seen in the [OIII]/H/3 map presented by'Ganda et al 1/200?), 
which shows a structure resembling the spiral pattern. Estimating 
ages and metallicity by comparison with SSP models, we infer 
that ages are mainly young-intermediate and the central metallicity 
is among the highest in our sample. The age and abundance ratio 
maps display the spiral arms structure as well, like the index maps. 
Fig. l21al (third column) demonstrates that for this galaxy constant 
star formation can be ruled out at a high confidence level: high r 
values are never likely. The r map shown in Fig.|6b](for age fixed to 
10 Gyr in the fitting procedure) presents a spatial structure that par- 
tially resembles the H/9 map and the spiral arms pattern. NGC 772 
is part of the sample studied in the series of papers by Cid Fer- 
nandes et al. 1 Cid-Fernandes et ^l2004l iGonzalez Delgado et al.l 
l2004ICid-Fernandes et al.ll2005h . who are engaged in a project to 
examine the stellar populations of low-luminosity active galactic 
nuclei. They classify this object as 'young-TO', i.e. a transition 
object between a LINER and an HII nucleus, with weak [O I] 
A6300 emission; they find that the nuclear spectrum is dominated 
by 10^-10^ yr-old stellar populations, with a contribution up to ~ 
32% from populations younger than 10^ yr and a similar contri- 
bution also from older stellar populations (^ 10® yr). Their data 
probe distances up to several hundred parsecs from the nucleus; 
they measure in NGC 772 significant radial variations of the stellar 



populations, with absorption line-strengths increasing in the inner 
3" indicating a younger population in the very centre, especially 
referring to the Ca II K line; the Mg I absorption strength presents 
instead less prominent variations over the spatial range of their 
data. Our data cover a different spectral range and a different 
spatial extension, but we also concluded that the populations in 
the inner parts have intermediate ages, and that the age (from the 
one-SSP analysis) increases in the inner ~ 5" (see Fig. lAlal ). 



NGC 864 

This is a barred Sc galaxy where a nuclea r radio source with th e 
linear size of » 300 p c has been detected jUlvestad & Holl2002h ; 
according to lEspada et al. (2005) its HI profile is symmetric in 
velocity but asymmetric in intensity; they also measure large-scale 
asymmetries in the two-dimensional kinematics, such as a warp. 
The mentioned asymmetries occur in the outer parts, well beyond 
the extent of the SAURON field. The asymmetries could have been 
induced by a past encounter with a companion, passing outside 
the optical disc but within the HI disc; the companion would have 
dispersed its stars and dark matter after the merging, consistent 
with the fact that NGC 864 appears to be remarkably isolated. 

Our line-strength maps show a very sharp central dip in both 
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Figure Alb. Same as Figure lATal for NGC 2964, NGC 3346, NGC 3423, NGC 3949, NGC 4030 and NGC 4 102. 



Fe5015 and Mgfc; in the region surrounding the dip, the Fe5015 
and Mgfc values rise. The H/3 absorption also takes on low values 
over an extended central region; as seen also in Table |2] the 
indices measured on the central aperture of this galaxy are among 
the lowest in the sample. From the one-SSP analysis, given the 
structure of the line-strength maps, we derive very low metallicity 
and old (~ 11. 22 Gyr) ages for the central aperture. 

The analysis of the /^y^ contours, in the continuous star 
formation approach, shows that this galaxy is unlikely to have 
experienced a constant star formation over time (see Fig. I21al 
fourth column). 



NGC 1042 

NGC 1042 is an Scd galaxy forming a pair with NGC 1035 at a 
separation of 22' (corresponding to 177 kpc). It has a bright, small 
nucleus and otherwise low surface brightness. Our line-strength 
maps are quite patchy, due to poor S/N in the data and consequent 
heavy binning of the spectra. The Mgfo index takes on low values 
in the very centre, and higher ones around it. H/3 is high in the 
centre and low outside it, being flat over most of the map; as 
can be seen also from the line-strength profiles in Fig. lAlal the 
line-strength indices stay approximately constant over most of the 
field, outside a central zone of radius ~ 3". The age of the stellar 
populations (one-SSP analysis) reaches a minimum in the centre. 



The metallicity is everywhere low, and peaks in the centre. 

The analysis of the Ax^ contours for the central aperture, in 
the continuous star formation approach, shows that this galaxy is 
not well described by an SSP (see Fig. I21al fifth column). The r 
map for age fixed to 10 Gyr (see Fig.l6c]l shows a peak in the very 
centre, corresponding to the peak in the H/3 and to the minimum in 
the yigb maps. 



NGC 2805 

NGC 2805 is an Sd galaxy seen nearly face-on and it is the brightest 
member of a multiple interacting system containing also NGC 2814 
(Sb), NGC 2820 (Sc) at 13' and IC 2458 (10). HI has been detected 
( lReakesllT97l iBosma et alj|l980t) and there are claims that the 
outer HI layers are warped (see for example |Bosma et al. 198(|). 
The galaxy seems to be also optically disturbed, since the spiral 
arms appear to be broken up into straight segments. 

The line-strength maps are rather patchy, displaying though 
a central peak in all three lines. By comparison of our measured 
indices with SSP models, we can state that the ages are young- 
intermediate over most of the field; that the metallicity peaks in the 
centre and that abundance ratios are super-solar. 

The analysis of the Ax^ contours, in the continuous star 
formation approach, shows that this galaxy is not very well 
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Figure Ale. Same as Figure lAlal for NGC 4254, NGC 4487, NGC 4775, NGC 5585, NGC 5668 and NGC 5678. 



described by an SSP, and that even a constant star formation rate 
cannot be excluded at all ages (see Fig. l21al sixth column). 



NGC 2964 

NGC 2964, classified as Sbc, forms a non-int eracting pair with th e 
10 galaxy NGC 2968 at 5! 8. As discussed by ' Canda et al] ( l2006h . 
the SAURON data reveal that this barred Sbc galaxy hosts an active 
galactic nucleus. When we calculate the line- strengths, we find 
that the Mgb index assumes low values in an elongated central 
region, higher values around it, and low values again towards 
the edges of the field; the Fe5015 index presents a very similar 
behaviour; the H/3 map has instead an opposite appearance, with 
an extended inner region with low index values and a surrounding 
area with higher values; in the very centre a shallow peak might be 
recognized. These features are recognizable also in the profiles in 
Fig. lAlbl We notice also that the central region with low Fe5015 
and Mgb corresponds to the ce ntral depression in t he emission-line 
[01Il]/H/3 map (see Fig. 5 in lOanda et Zll2006h . possibly indi- 
cating ongoing star formation. Converting the line-strength maps 
to age and metallicity via our comparison with the SSP models, 
we find young-intermediate ages over most of the field, with a 
peak in the central region, and low metallicity everywhere, with a 
minimum in the central region, a positive gradient out to « 5", and 
then a rather constant behaviour out to ~ 15". From the continuous 



star formation approach, we cannot exclude at all ages the case 
of constant star formation, but the r map for age fixed to 10 Gyr 
(in the fitting procedure) presents a large central region with low 
values, possibly suggesting a « 5 Gyr long starburst. 



NGC 3346 

This Scd galaxy is very poorly studied in the literature. Ages of 
the stellar populations (one-SSP approach) are young-intermediate 
and metallicity low over most of the field. From the continuous 
star formation approach, we cannot exclude at all ages the case of 
constant star formation (see Fig. l21bl second column). 



NGC 3423 

This Sbc galaxy is also poorly known. Both the Fe5015 and Mgb 
indices peak in an extended central region and decrease towards the 
edge of the field; the H/3 map presents opposite behaviour. The r 
map (for age fixed to 10 Gyr in the fitting procedure) displays low 
values over a large central area, possibly indicating a region where 
the bulk of star formation happened in a quite short starburst. 
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NGC3949 

NGC3949 is an Sc galaxy belonging to the Ursa Major cluster. 
NGC3949 presents quite regular Fe5015 and Mgb line-strength 
maps, peaking in an extended central region and decreasing rather 
smoothly moving outwards. The H/3 index is low in the central 
region and higher in the surroundings; interestingly, it reaches a 
minimum in a zone south of the centre where the emission-line 
ratio [0III]/H/3 is particularly depressed (see relative Figure in 
iGanda et al.l |20og). The values measured for the H/3 index are the 
highest in the sample, as one can see from Table [2] and from Fig. 
lAlbl Ages (from the one-SSP approach) are young (~ 1.0 Gyr) in 
the central region, somewhat higher around it and decreasing again 
towards the edge of the field. The metallicity is everywhere low, 
and decreasing outwards. 

Constant star formation cannot be ruled out; it is likely to 
characterise the star formation mode over most of the field (see 
the r map for age fixed to 10 Gyr in Fig. [6el r assumes high 
values). Actually, when performing the fit allowing the code to 
choose among exponentially declining and constant star formation, 
constant is selected over most of the field. 



NGC4030 

This Sbc galaxy is characterised by v ery regular kinemat ics for both 
the stellar and gaseous components dGanda et alJI 20061) . The line- 
strength maps also have a smooth appearance, with Fe5015 and 
Mgb decreasing at larger radius and H/3 increasing. The H/3 map 
seems to be asymmetric with respect to the centre, since the western 
side displays higher values than the eastern one. The stellar popula- 
tions are young over the whole field and the metallicity, among the 
highest in our sample (the third after NGC 488 and 772), decreases 
moving outwards; the abundance ratio is approximatively solar. 
The r map for age fixed to 10 Gyr in Fig. [6fl shows a large central 
depression. The analysis of the Ax^ contours, in the continuous 
star formation approach, shows that this galaxy is probably not very 
well described by an SSP, nor by a constant star formation rate. 

NGC 4102 

NGC 4102 is an Sb galaxy, classified as a LINER in the NED 
database; it is kn own to be a powerful far-infrared galaxy 
I Young et a l. 1989) and also to h ave a strong nuclear radio source 
I Condon et al. 1982 ). Devereuxl {l989) classified it as one of the 
most powerful nearby starburst galaxies. The sign s of activity 
are d etectable also in the SAURON data: as noted bv lGanda et al.l 
( I200 6). the [OIII] maps trace the outflowing gas, with a region of 
very high gas velocity dispersion corresponding also to high values 
in the [0III]/H/3 line ratio map, as expected for an active object. 
In that same region, we observe double-peaked line profiles for the 
[OIII] emission-lines, pointing towards the presence of different 
components in the emitting material. As for the line-strength 
indices, Mgb is low in an inner region, increases around it, and 
decreases again in the outer parts. The Fe5015 map is instead more 
centrally concentrated. The map of H/3 absorption looks more 
complex and not symmetric. Ages (from our one-SSP analysis) 
are young all over the field and metallicity is low. Neither in the 
line-strength nor in the age and metallicity map is there a clear 
connection with the mentioned region of activity. On the basis of 
our continuous star formation analysis, this galaxy is probably not 
very well described by an SSP. 



NGC 4254 

Also known as M99, NGC 4254 is a bright Sc galaxy located on 
the periphery of the Virgo cluster, at a projected distance of 3.7° 
(~ I Mpc) from the centre of the cluster. Its optical appearance is 
dominated by a peculiar one-arm structure: the arms to the north- 
west are much less defined than the southern arm. This kind of 
spiral structure could be related to an external driving mechanism, 
but for NGC 4254 there are no close companions. Phookun, Vogel 
& Mundy (1993) carried out deep HI observations of NGC 4254, 
detecting non-disc HI clouds, with velocities not following the disc 
velocity pattern and an extended low surface density tail northwest 
of the galaxy; they interpreted these observational results as 
infall of a disintegrating gas cloud. They also noticed that there 
is an hole in the HI emission at the centre of the galaxy, with a 
diameter of ~ 3 kpc. VoUmer, Huchtmeier & van Driel (2005) 
propose instead a scenario where NGC 4254 had a close and rapid 
encounter with another massive galaxy when entering the Virgo 
cluster: the tidal interaction caused the one-arm structure and the 
HI distribution and kinematics is due to ram pressure stripping. 
Recently, , Soria & Wong 2006, studied the X-ray properties of M99, 
noticing that a phenomenon often associated with tidal interactions 
and active star formation -both documented- is the presence of 
ultraluminous X-ray sources (ULXs). The X-ray emission appears 
approximatively uniformly diffused across the inner ~ 5 kpc from 
the nucleus. They do find a ULX, amongst the brightest observed, 
located at ~ 8 kpc southeast of the nucleus, close to the position 
where a large HI cloud seems to join the gas disc. From observa- 
tions of the radial velocities, they suggest that the cloud is falling 
onto the galactic disc and try to build a link between this collisional 
event and the formation of the bright ULX. Our data cover a much 
smaller spatial extent than the mentioned observations, therefore a 
comparison is not possible; within our field, both the Fe50I5 and 
Mgb maps show a regular structure with a central concentration; 
at the northern edge of the field, both Fe50I5 and Mgb seem to 
increase again; the H/3 map has opposite behaviour, presenting a 
central depression. The ages (from our one-SSP approach) are gen- 
erally young and the metallicity, quite low everywhere, decreases 
slightly moving outwards and shows the same enhancement at the 
northern edge of the field that we noticed for the Fe50I5 and Mgb 
line-strength. The r map (see Fig. \6g\ shows a central depression. 
The analysis of the Ax^ contours, in the continuous star forma- 
tion approach, shows that NGC 4254 is probably not very well 
described by an SSP, nor by a constant star formation rate. 

NGC 4487 

This Scd galaxy forms a pair with the Scd galaxy NGC 4504, at a 
separation of 35' (corresponding to ~ 165 kpc). The Fe5015 and 
Mgb maps seem to be low in the very centre and higher outside 
it; from the profiles in Fig. lAIcl one can see that outside the inner 
circle of ~ 3" radius, the line-strengths keep an approximately 
constant value, out to ~ 13". Ages (one-SSP approach) are 
intermediate-young and metallicities very low, with respectively a 
maximum and a minimum in the centre. Constant star formation 
cannot be ruled out. 

NGC 4775 

NGC 4775 is an Sd galaxy, for which very little is known from 
previous studies. The metal lines (Fe50I5 and Mgb] seem to be 
centrally concentrated, while H/3 is low over an extended central 
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area, and higher in the surrounding region. The maps are not 
synmietric with respect to the centre. The galaxy is young in the 
nucleus and has older ages in an annular (extended) region around 
the centre; the metal content is everywhere low, and peaks in the 
centre. Constant star formation cannot be ruled out, and is selected 
over a large part of the field when allowing to choose among 
constant and exponentially declining star formation. 

NGC5585 

NGC5585 is a barred Sd galaxy; together with NGC5204 (Sm), 
NGC5474 (Scd), NGC5477 (Sm) , HoIV (Im) and MlOl (Scd) it 
forms the MlOl group. The Fe5015 and Mgb indices assume val- 
ues among the lowest in the sample; H/3 has instead a value among 
the highest ones (see Table |2}. The metallicity retrieved from the 
one-SSP analysis is everywhere very low. In the continuous star 
formation approach, constant star formation cannot be ruled out; in- 
deed, when allowing to choose among exponentially declining and 
constant star formation, in the minimisation, for age fixed to 10 
Gyr the constant star formation is selected (for the central aperture). 

NGC5668 

NGC 5668 is another Sd galaxy, which has a high rate of star for- 
mation, as indicated by its large far-infrared and Hq luminosities 
jSchulman et alj [l996l and references therein). The metal lines 
(Fe5015 and Mgb) appear to be centrally concentrated, while H/3 
has a more complex appearance. The stellar populations are young 
and metal poor all over the field; the age rises out to ~ 8" and 
then decreases (see Fig.[6i]for a full view of the spatial structure of 
the age map and Fig. lAlcl for the azimuthally averaged profiles). 
Constant star formation cannot be ruled out at all ages, and is 
selected over a large part of the field when allowing to choose 
among constant and exponentially declining star formation. 

NGC 5678 

This is a barred and very dusty Sb galaxy. Fe5015 and Mgb 
have a quite regular and centrally concentrated appearance; H/3 
is depressed in the inner region and increases in a annulus-like 
structure around it. In the very centre H/3 has a local maximum, 
as shown also in the azimuthally-averaged profile in Fig. lAlcl 
Ages (one-SSP approach) are young over most of the field and the 
metallicity, rather low everywhere, decreases moving outwards. 
The T map at age fixed to 10 Gyr (see Fig. |6i} shows a large 
central depression. The analysis of the Ax^ contours shows that 
NGC 5678 is probably not very well described by an SSP, nor 
by a constant star formation rate. Like NGC 628 and 772, this 
galaxy is part of the sample stud i ed in the Cid Femandes et al 
papers (Cid-Fernandes et al. l l2004lGonzaiez Delgado et aLl l2004 
[Cid -Fernandes et al. 200|); following their classification, it is 
a 'young-TO' galaxy, with weak emission in the [O I] A6300 
line; according to those authors, the nuclear spectrum is heavily 
dominated by 10**-10^ yr old stellar populations, with negligible 
contribution from stars younger than lO'^ yr. They detect radial 
variations of the absorption indices, particularly in the Ca II K line, 
but less clearly evident than in the case of NGC 772. 



